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THE  TURBULENT  BOUNDARY  LAYER  ON  A  CYLINDER  IN  AXIAL  FLOW 

1.  INTRODUCTION 

The  effect  of  transverse  curvature  on  the  turbulent  boundary  layer  that  develops  as  a  fluid  flows 
parallel  to  a  cylindrical  surface  has  applications  in  several  different  fields.  Early  interest  in  the  effect  of 
tiansverse  curvature  resulted  from  efforts  to  develop  frictional  resistance  similarity  laws  to  allow  scaling  from 
mode  s  to  full-scale  for  ship  hulls.1  Other  early  interest  was  in  the  boundary  layer  that  develops  on  cylindrical 
me  fes2  or  during  the  towing  of  submerged  cables  or  cylindrical  bodies.3  Finally,  interest  in  the  drag  of 
glas  cr  polymer  fibers  moving  through  air  during  fabrication  motivated  early  research.4 

Like  the  planar  boundary  layer,  the  cylindrical  boundary  layer  is  two-dimensional  in  the  streamwise 
a;  wi-i-normal  directions,  while  it  is  periodic  in  the  spanwist  direction.  However,  most  boundary  layer 
resea.  an  has  been  concentrated  on  the  planar  boundary  layer.  In  spite  of  the  apparent  simplicity  of  an 
a-isyn  ,ieiric  boundary  layer,  the  volume  of  research  and,  hence,  the  level  of  understanding  regarding  the 
pi  ysics  O’  trns  boundary  layer  have  lagged  behind  that  of  planar  boundary  layers. 

Jpcr  further  investigation,  an  axisymmetric  boundary  layer  is  not  as  simple  as  its  two-dimensional 
tr.a;  icier  implies  because  cf  the  existence  of  an  additional  length  scale,  the  radius  of  transverse  curvature, 
a.  ~\  nis  aoditiona  length  scale  suggests  several  nondimensional  parameters,  i.e.,  a+=  a  U^v ,  5/a  ,  Rg=  a 
L',fc.-  ,  a-  ;  ^=W(  v  1  U^a2).  The  first  scaling  parameter,  a+,  is  a  wall  scaling  based  on  the  friction  velocity,  U^, 
amx,  the  inematic  v  scosity,  v.*  The  second  scaling  parameter  is  a  transverse  curvature  ratio  relating  the 
sc  unda  y  layer  tlx  xness,  8,  with  the  radius  of  transverse  curvature  of  the  surface.  The  third  scaling 
p  ara  neter  mixes  trn.  outer  free  stream  velocity,  U„,  with  the  transverse  curvature  resulting  a  Reynolds 
number  based  on  transverse  curvature,  Ra.  The  last  scaling  parameter  incorporating  the  streamwise 
coordinate,  x,  was  proposed  in  conjunction^  with  the  laminar  axisymmetric  boundary  layer.6  Note  that  all  of 
these  scaling  parame  a  except  Ra  are  "local"  with  respect  to  the  flow.  In  other  words,  they  depend  upon 
Ux  ,  8,  or  x  which,  In  tu~  depend  upon  the  streamwise  position  in  the  boundary  layer.  Only  Ra  is 
independent  of  tt  sirearriwise  location  in  the  boundary  layer. 

*  All  va  tables  subscripted  with  +  have  been  nondimensionalized  using  v  and  Ux=V(xw/p), 
where  xw  is  the  w;  .  she  or  stress  and  p  is  the  fluid  density. 

t  Seban  ar  j  Bond6  indicate  that  the  parameter  \  was  originally  proposed  by  Young,6  although 
this  reference  was  riot  available  to  the  author  for  confirmation. 


1 


TR  8389 


Consider  the  extreme  cases  of  transverse  curvature.  Firot,  as  the  radius  of  transverse  curvature, 
a,  approaches  infinity  (a+  and  Ra  very  large;  5/a  and  %  very  small),  the  boundary  layer  should  be  similar  to  a 
planar  boundary  layer.  A  similar  result  occurs  when  the  boundary  layer  thickness  is  very  small  compared  with 
the  radius  of  transverse  curvature  resulting  in  small  8/a.  On  the  other  hand,  if  the  radius  of  transverse 
curvature,  a,  is  very  small  (  a+  and  Ra  small;  5/a  ar.d  \  large)  the  boundary  layer  can  be  much  thicker  than  the 
diameter  of  the  cylinder.  For  instance,  Richmond7  measured  the  cylindrical  boundary  layer  to  be  nearly  an 
inch  (2.5  cm)  thick  on  a  wire  with  a  radius  of  0.012  inches  (0.03  cm).  Denli  and  Landweber8  note  that  as  5/a 
increases,  the  outer  flow  in  a  cylindrical  boundary  layer  becomes  increasingly  independent  of  the  wall 
suggesting  a  flow  that  would  be  similar  to  a  cylindrical  wake  with  a  modified  inner  boundary  condition. 

The  range  of  experimental  data  that  is  available  for  the  mean  velocity  profile  is  shown  in  figure  i  in 
of  terms  various  scaling  parameters,  in  addition,  the  Reynolds  number  based  on  momentum  thickness,  R0, 
is  included,  since  this  parameter  is  also  defined  for  planar  boundary  layers.  The  range  of  measurements  for 
a+  and  Ra  span  three  orders  of  magnitude,  but  the  same  measurements  only  represent  a  range  of  one  or  two 
orders  of  magnitude  for  5/a  and  £,.  Luxton  et  al.9  show  a  logarithmic  correspondence  between  a+  and  Ra 
for  available  experimental  data.  *  Their  figure  also  shows  that  as  a+  and  Ra  increase,  5/a  decreases  for  the 
available  expenmental  data.  The  direct  relation  between  a+  and  Ra  is  primarily  a  result  changing  the  radius 
Of  transverse  curvature  or  the  free  stream  velocity  to  obtain  a  variation  in  scaiing  parameters.  Usiny  i’nis 
approach,  a+,  5/a,  and  Ra  are  not  independently  varied.  Very  little  data  are  available  based  on  parameter 
variation  accomplished  by  measuring  the  boundary  layer  characteristics  at  different  axial  locations  given  the 
same  transverse  curvature  in  order  to  vary  a+  and  5/a  while  keeping  Ra  constant.  Because  data  obtained  in 
this  way  are  scarce,  proper  scaling  parameters  have  been  difficult  to  identify. 

At  this  point  it  is  necessary  to  discuss  a  nuance  in  the  terminology  that  is  often  overlooked.  The 
terms  "axisymmetric  boundary  layer”  and  "cylindrical  boundary  layer"  have  been  used  interchangeably  to 
describe  a  boundary  layer  that  develops  as  a  fluid  (lows  parallel  to  the  axis  ot  a  cylinder.  As  will  become 
evident  shortly,  obtaining  axial  symmetry  experimentally  is  very  difficult.  The  problem  becomes  increasingly 
acute  as  the  boundary  layer  becomes  thick  with  respect  to  the  radius  of  the  cylinder.  Because  the  boundary 
layer  that  develops  on  a  cylinder  in  axial  llow  is  not  necessarily  axisymmetric,  the  term  "cylindrical"  instead  of 
"axisymmetric’'  will  be  used  to  describe  this  boundary  layer. 

"The  labels  on  the  axes  of  figurel  of  Luxton  et  al.9  are  inadvertently  reversed  and  should  be 
interchanged. 
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Fig.  1  Range  ol  Values  of  Nondimensional  Parameters  for  Which  Mean  Velocity  Data 
for  a  Cylindrical  Boundary  Layer  Are  Availaole 
(LBR  =  Luxton  et  al.  9  LLS  =  Lueptow  et  al.,19  RK  =  Rao  and  Keshavan,11 
R=Richmond/  WWSB  =  Willmarth  et  al.12^ 
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2.  MEAN  VELOCITY  PROFILE 

The  differences  between  the  mean  velocity  profiles  for  a  cylindrical  boundary  layer  and  a 
planar  boundary  layer  are  substantial,  especially  as  the  radius  of  transverse  curvature  becomes  small. 
This  is  evident  in  the  velocity  defect  data  of  Willmarth  et  al.1 2  reproduced  in  figure  2.  For  the  larger 
cylinders  (corresponding  to  a  small  5/a )  the  velocity  defect  profile  is  nearly  identical  to  the  profile  for  a  flat 
plate.  As  the  cylinder  diameter  becomes  smaller  (corresponding  to  larger  5/a ),  the  velocity  defect 
profiles  become  fuller.  Richmond7  and  Willmarth  and  Yang"14  suggest  that  this  characteristic  is  similar  to 
that  for  a  laminar  axisymmetric  boundary  layer.  Glauert  and  Lighthill15  note  that  in  a  laminar  cylindrical 
boundary  layer  the  shear  force  per  unit  length  on  a  cylinder  equals  the  circumference  times  the  shear 
stress,  and  that  this  force  is  the  same  on  all  cylindrical  fluid  surfaces  near  the  wall  (see  equation  2-1). 
Consequently,  the  shear  stress  must  vary  with  1/r,  where  r  is  the  radial  coordinate.  Since  the  shear 
stress  is  proportional  to  5u/3r  (where  U  is  the  mean  velocity  in  the  streamwise  direction),  then  0U/3r 
must  also  vary  with  1/r.  As  the  cylinder  radius  becomes  smaller,  1/r  can  become  larger  so  the  velocity 
gradient  increases.  Thus,  a  fuller  velocity  profile  results  as  the  radius  of  the  cylinder  decreases. 


The  differences  between  a  planar  boundary  layer  and  a  cylindrical  boundary  layer  are  also 


evident  when  the  mean  velocity  is  plotted  as  a  function  of  the  distance  from  the  wall  in  the  usual  wall 


coordinates,  U+  and  y+  as  shown  in  figure  3.  As  5/a  increases,  the  cylindrical  mean  velocity  profiles 
drop  below  the  planar  profile.  For  small  5/a,  the  log  region  of  the  mean  velocity  profile  nearly  mate*-  -  - 


the  log  region  for  a  pianar  boundary  layer. 


2.1  MEASUREMENTS  OF  THE  WALL  SHEAR  STRESS 


Although  it  may  seem  premature  to  describe  measurements  of  the  wall  shear  stress  at  this 
point,  it  is  necessary  to  do  so  because  the  friction  v'/ocity  appears  in  the  scaling  of  the  mean  velocity 
profiles  when  presenting  the  mean  velocity  profile  in  wall  coordinates.  However,  a  major  weakness  in 
the  presentation,  analysis,  and  use  of  nearly  all  of  the  velocity  data  for  cylindrical  boundary  layers  is  the 
difficulty  in  measuring  the  wall  shear  stress,  .  Although  techniques  have  been  developed  for  use  in 
measuring  the  wall  shear  stress  in  a  planar  boundary  layer,  these  methods  are  not  directly  applicable  to  a 
cylindrical  boundary  layer  because  of  the  effect  of  transverse  curvature. 

Richmond7  made  an  attempt  to  measure  the  skin  friction  directly  on  a  0.25  inch  (0.64  cm) 
diameter  cylinder  using  a  floating  element  with  an  inductance  coil  position  transducer.  The  measured 
coefficient  of  friction  for  hypersonic  flow  is  higher  than  that  for  a  planar  boundary  layer  at  the  same 
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momentum  thickness  Reynolds  number.  The  skin  friction  was  not  measured  for  incompressible  flow. 

In  several  studies  the  coefficient  of  drag  of  a  cylinder  in  axial  flow  was  measured  directly  These 
experimental  studies  fall  into  three  categories.  The  first  category  is  that  of  towing  experiments  where  a 
cylindiical  body  of  finite  length  is  towed  through  a  fluid.  Unfortunately,  in  several  of  these  tests  it  is 
unclear  whether  the  flow  was  laminar  or  turbulent.  Kempf 1 7  and  Hughes1 8  towed  neutrally  buoyant 
cylinders  up  to  84  ft  (26  m)  long  in  a  water  tow  tank  and  measured  the  total  drag.  These  results  showed 
that  the  coefficient  of  friction  decreases  as  a  function  of  both  the  length-to-diameter  ratio  and  the  length 
Reynolds  number,  Rl-U^L/v,  where  l  is  the  length  of  the  cylinder.  Zajac1®  carried  out  similar  tests  on 
smooth-  and  rough-surfaced  cables  that  were  not  neutrally  buoyant,  although  the  length-to-diameter 
ratio  was  not  varied  (see  reference  3  for  these  results  in  terms  of  the  diameter  Reynolds  number). 
Towing  tests  of  optical  fibers  and  monofilaments  were  carried  out  by  Kennedy  et  al.20  By  using  fibers  as 
small  as  0.2  mm  diameter  and  60  to  70  m  long,  the  drag  coefficient  for  very  large  length-to-diameter 
ratios  and  a  wide  range  of  diameter  Reynolds  numbers  were  measured.  Their  results  indicate  that  the 
average  coefficient  of  friction  decreases  as  the  diameter  Reynolds  number  or  the  length-to-diameter 
ratio  increases.  However,  because  of  the  very  small  diameter  of  the  fibers  that  were  used,  the  wake  of 
the  tow  body  and  support  stmt  may  have  influenced  the  drag  over  the  portion  of  the  fiber  immediately 
behind  them: 

In  the  second  category  of  direct  drag  measurement,  cylinders  were  suspended  in  a  wind  tunnel 
and  the  coefficient  of  drag  was  measured  directly.  Again  the  degree  to  which  the  flow  field  was 
turbulent  in  these  tests  is  difficult  to  ascertain.  Andrews  and  Cansfield21  measured  the  drag  force  per 
unit  length  in  this  way  although  experimental  details  are  sketchy.  Gould  and  Smith22  measured  the 
coefficient  of  drag  on  monofilaments  ranging  from  8.5  p.m  to  150  )i.m  in  diameter  and  91  cm  long.  At 
higher  diameter  Reynolds  numbers  the  length-to-diameter  ratio  had  little  effect  on  the  drag  coefficient. 
However,  ax  lower  diameter  Reynolds  numbers,  the  drag  coefficient  appears  to  have  some  dependence 
upon  the  length-to-diameter  ratio,  Ni  and  Hansen23  estimated  the  drag  on  a  flexible  cylinder  in  a  water 
tunnel  by  measuring  the  elongation  of  the  cylinder  and  correlating  that  elongation  to  an  equivalent  axial 
load. 


The  third  category  of  direct  drag  measurement  has  been  in  connection  with  the  tension  of 
polymer  or  glass  fibers  during  production.  Typically,  the  tension  of  a  continuously  extruding  fiber  is 
measured  using  a  tensiometer.  The  aerodynamic  drag  on  a  fiber  is  found  by  subtracting  the  effects  of 
gravity,  rheological  drag,  surface  tension,  and  inertia  from  the  total  tension  force.24  Because  the  fibers 
are  very  small,  the  transverse  curvature  is  certain  to  play  a  large  role.  In  some  cases,  6/a  can  be 
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estimated  to  be  as  large  as  of  the  order  of  200. 25  The  results  of  several  experiments  of  this  type  show 
that  the  drag  coefficient  decreases  with  increasing  radius  Reynolds  number,  Rea26"29  Unfortunately, 
these  drag  coefficient  measurements  are  subject  to  question  for  two  reasons.  First,  it  is  often  unclear 
whether  the  flow  is  laminar,  transitional,  or  turbulent.  Second,  the  effect  of  the  transverse  vibration  of 
the  filament  is  unknown. 

An  indirect  method  of  evaluating  the  wall  shear  stress  in  cylindrical  boundary  layers  has  been 
through  the  use  of  a  Preston  tube,  a  circular  impact  pressure  tube  in  contact  with  the  wall.  Willmarth  et 
al.12  found  that  as  long  as  the  cylindrical  boundary  layer  profile  coincides  with  the  planar  profile  the  wall 
shear  stress  can  be  measured  using  a  Preston  tube,  although  the  Preston  tube  diameter  must  be  less 
than  0.3  times  the  cylinder  diameter.  Using  this  method,  they  found  fhat  the  coefficient  of  friction 
generally  decreases  with  increasing  Rea,  Re0,  and  a+,  but  is  consistently  larger  than  the  coefficient  of 
friction  for  a  flat  plate  boundary-  layer.  Similar  results  were  reported  by  Willmarth  and  Yang14  and  Yu.30 
For  very  small  cylinders  Willmarth  et  al.  were  not  confident  of  results  obtained  using  a  Preston  tube 
because  of  differences  between  the  planar  velocity  profile  and  the  cylindrical  velocity  profile. 

The  inherent  weakness  in  the  use  of  a  Preston  tube  is  that  the  calibration  data  for  a  flat  plate 
must  be  assumed  valid  for  a  cylindrical  boundary  layer.  Since  a  Preston  lube  extends  well  into  the  log 
region  of  the  boundary  layer,  and  the  log  region  of  a  cylindrical  boundary  layer  is  significantly  different 
from  that  of  a  planar  boundary  layer  for  large  5/a,  the  validity  of  the  flat  plate  calibration  is  questionable  for 
small  diameter  cylinders.  To  avoid  this  problem,  Lueptow  and  Haritonidis31  used  a  modified  impact 
pressure  p:obe  so  small  that  the  opening  was  within  the  sublayer.  Because  of  the  small  difference  in  the 
sublayer  velocity  profiles  for  the  planar  case  and  the  cylindrical  case,  they  assumed  that  the  planar 
calibration  could  be  used  for  the  cylindrical  boundary  layer  measurements.  Using  this  method,  they 
tourid  the  wall  shear  stress  to  be  greater  than  the  pianar  boundary  iayer  values  and  to  decrease  with 
increasing  Rex. 

Several  indirect  measures  of  the  wall  shear  stress  have  been  used.  The  most  common  method 
is  the  use  of  a  Clauser  plot13  to  match  the  slope  of  the  log  region  of  the  mean  velocity  profile  to  that  for 
the  planar  case  (see  for  example  reference  32).  This  method  will  certainly  provide  inappropriate  results 
when  the  radius  transverse  curvature  is  small,  because  the  planar  and  cylindrical  log  region  profiles  are 
known  to  be  different. 

Willmarth  et  al 12  and  Lueptow  et  al.10  used  a  method  in  which  the  measured  mean  velocity 
|  profile  data  in  the  sublayer  and  buffer  zone  were  fitted  to  the  planar  sublayer  profile  through  the 
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adjustment  of  the  friction  velocity,  thus  specifying  the  wall  shear  stress.  The  assumption  here  was  that 
very  c  lose  to  the  wall  the  differences  between  a  planar  profile  and  a  cylindrical  profile  are  negligible. 
Willm,  rth  et  al.  verified  this  assumption  for  5/a  <10  by  comparing  the  wall  shear  stress  obtained  using 
this  method  with  the  wall  shear  stress  measured  using  a  Preston  tube.  However,  Lueptow  and 
Haritonidis31  found  that  matching  the  mean  velocity  profile  to  the  sublayer  profile  resulted  in  wall  shear 
stresses  that  were  higher  than  those  measured  with  their  sublayer  impart  pressure  probe. 

Lueptow  et  a!.,10  noting  the  relationship  between  the  Reynolds  stress  and  the  wall  shear 
stress  in  a  cylindrical  boundary  layer  (see  equations  2-1  and  2-2),  used  the  slope  of  the  Reynolds  stress 
profile  when  plotted  as  a  function  of  a/r  to  estimate  the  wall  shear  stress.  Although  this  method  shows 
trends  consistent  with  other  methods,  the  scatter  in  the  data  is  too  large  to  make  this  a  useful  technique 
for  finding  the  wall  shear  stress. 

Finally,  it  ;hould  be  noted  that  the  wall  shear  stress  can  be  found  from  the  streamwise  gradient 
of  the  momentum  integral  expression  for  an  axisymmet.ic  boundary  layer  after  substitution  of 
experimental  velocity  profiles  into  the  integral  (see  equation  2-7).  Afzal  and  Singh32  used  this  method 
and  obtained  resuhs  similar  to  ttnse  obtained  using  Clauser  plots  for  small  5/a.  In  practice,  this 
procedure  is  erroneous  unless  perfect  axisymmetry  is  maintained  as  pointed  out  by  Lueptow  et  al/*3 
Errors  related  to  the  differentiation  of  experimentally  derived  quantities  is  also  a  problem  with  this 
method. 


2.2  COMMENTS  ON  EXPERIMENTAL  METHODS  FOR  VELOCITY 
MEASUREMENTS  IN  CYLINDRICAL  BOUNDARY  LAYERS 

Most  experimental  setups  used  to  make  velocity  measurements  of  cylindrical  boundary  layers 
involve  the  suspension  of  a  cy'!ndrical  model  along  the  centerline  of  a  wind  tunnel  and  the  measurement 
of  the  velocities  in  the  boundary  layer  using  hot  <vire  anemometry.  This  approach  brings  about  a  host  of 
exDerirnental  problems.  Most  crucial  perhaps,  is  the  problem  of  maintaining  axisymmetry  of  the 
boundary  layer.  In  spite  of  holding  i he  cylinder  in  tension9  or  supporting  it  with  guy  wires,30  asymmetry 
of  the  boundary  layer  often  results  from  the  sag  of  the  cylinder  in  a  horizontal  wind  tunnel.  The  sag 
problem  has  been  successfully  overcome  by  the  use  of  a  vertical  wind  tunnel  by  Willmarth  et  al.12  Even 
so,  they  experienced  difficulty  in  finding  cylinders  that  ware  straight.  Apart  trom  the  sag  problem  is  the 
di  ficulty  in  keeping  the  cylinder  parallel  to  the  air  fiow  in  the  wind  tunnel.  Willmarth,  Shama  and  Inglis33 
a'tf  Lueptow  et  al.10  found  that  an  angle  as  small  as  0.05r  between  the  axis  of  the  cylinder  and  the  f!ov 
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direction  can  alter  the  boundary  layer  thickness  substantially.  Rao  and  Keshavan1 1  investigated  the 
growth  of  the  momentum  thickness  with  streamwise  location  for  a  yaw  angle  of  1°. 


The  leading  edge  condition  is  also  troublesome  for  the  experimentalist.  Often  a  cylinder  is 
mounted  in  a  wind  tunnel  so  that  the  cylinder  extends  upstream  into  the  contraction  section,9  the 
settling  chamber,12  or  outside  of  the  wind  tunnel10  to  create  an  intentionally  ambiguous  leading  edge 
condition.  Of  course  this  leads  to  problems  in  detining  any  streamwise  parameters  (for  example  see 
reference  31 ).  Sometimes  it  is  necessary  to  trip  the  boundary  layer  in  order  to  obtain  turbulent  flow 
within  the  test  section.  Lueptowet  al.10  simply  used  an  O-ring  around  the  cylinder  at  the  upstream 
end  of  the  test  section  tor  this  purpose.  A  significant  number  of  experiments  have  been  performed  with 
a  cylinder  moving  through  a  quiescent  fluid  (for  example  reference  25).  Again  this  experimental  setup 
results  in  ambiguity  in  the  definition  of  streamwise  parameters  and  the  leading  edge  conditions.  !n  cases 
where  the  leading  edge  of  the  cylinder  is  within  the  test  section,  hemispherical30  and  ogive-shaped1 1 
leading  edge  geometries  have  been  used. 


Velocity  measurements  in  boundary  layers  on  very  small  cylinders  are  difficult  because  the 
length  scale  of  the  hot-wire  probe  or  pitot  tube  is  often  of  the  same  order  as  the  radius  of  curvature.  For 
instance,  in  widely  referenced  experiments,  Richmond7  used  a  0.015  inch  (0.38  mm)  long  hot-wiro  to 
make  mean  velocity  measurements  on  a  cylinder  with  a  radius  of  only  0.012  inch  (0.30  mm).  Near  the 
wall,  the  measured  velocity  was  probably  in  error  since  the  hot-wire  averages  the  velocity  over  the 
length  of  the  wire,  and  the  ends  of  the  hot-wire  were  exposed  to  higher  mean  velocities  than  the  center 
of  the  wire.  In  addition,  the  midpoint  of  the  trot-wire  must  be  centered  over  the  cylinder,  so  ‘hat  it  will  not 
be  exposed  to  different  velocities  at  each  end  of  the  wire. 


In  the  sublayer  region,  the  analysis  of  the  mean  velocity  profile  begins  with  the  streamwise 
Navier-Stokes  equation  in  cylindrical  coordinates.  Assuming  ihere  is  no  pressure  gradient  and  noting 
that  the  streamwise  gradients  and  the  wall-normal  velocity  are  small  in  the  sublayer,  the  equation 
reduces  to 


a\v  “  rx- 

The  total  shear  stress,  t,  is  given  by 


(2-1) 
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T/p  =  V  (3U/dr)  -  uv  , 


(2-2) 


where  p  is  the  fluid  density,  u  and  v  are  the  fluctuating  velocities  in  the  streamv/ise  and  wall-normal 
directions,  respectively,  and  the  overbar  denotes  time  averaging.  The  second  term  on  the  right-hand 
side  of  (2-2)  is  the  Reynolds  stress  which  arises  from  the  turbulent  velocity  fluctuations.  Equation  (2-1), 
sometimes  called  the  constant  shear  moment,  is  analogous  to  the  constant  shear  layer  very  near  the 
wall  in  a  planar  boundary  layer. 

Reid  and  Wilson3  and  Rao34  reasoned  that  the  fluctuating  velocities  are  small  in  the  viscous 
sublayer  so  the  Reynolds  stress  is  negligible  allowing  the  integration  of  (2-1)  after  substituting  in  (2-2). 
This  leads  to 

U+»a+ln(r/a).  (2-3) 

Equation  (2-3)  is  widely  accepted  as  the  appropriate  description  of  the  mean  velocity  profile  in  the 
sublayer  of  a  cylindrical  boundary  layer.  Noting  that  !ri(r/a)  =  ln[1  +(y/a)]  =>  y/a  tor  smali  y/a,  (2-3)  reduces 
to  the  planar  sublayer  expression,  U+  *  y+  ,  when  the  radius  of  curvature,  a,  is  large. 

Experimental  verification  of  equation  (2-3)  is  difficult  since,  like  the  planar  sublayer  expression, 
this  expression  is  valid  only  out  to  a  distance  from  the  wall  of  y+<0(10 ).  Velocity  measurement  is  very 
difficult  this  close  to  a  surface  because  of  heat  transfer  to  the  wall,  probe  resolution,  and  accurate 
calibration  at  very  low  velocities.  Nevertheless,  Willmarthet  al. 12  showed  that  very  near  the  wall  Jy+*5) 
the  data  appear  to  fall  on  the  curve  defined  by  equation  (2-3)  as  shown  in  figure  3.  In  that  same  figure,  it 
is  evident  that  equation  (2-3)  is  not  significantly  different  from  the  planar  sublayer  profile,  oven  tor  the 
smallest  a+  that  Wiiimarth  et  al.  measured. 

2.4  MEAN  VELOCITY  PROFILE  IN  THE  INNER  REGION 

Away  from  the  wall,  outside  the  viscous  sublayer,  is  the  inner  region  of  flow  where  the 
influence  of  the  wall  is  important,  yet  viscous  shear  stresses  do  not  dominate.  The  equivalent  region  in 
a  planar  boundary  layer  is  often  called  the  log  region  Many  attempts  have  been  made  at  the 
measurement  of  the  mean  velocity  profile  in  the  inner  region  as  indicated  in  table  1 .  Several  of  these 
studies  were  flawed  by  insufficient  probe  resolution  as  noted  in  the  table.  Of  course,  ihe  problem  of 
probe  resolution  is  exaggerated  for  the  very  situation  where  the  effects  of  transverse  curvature  are 
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greatest  -  for  small  cylinders  where  the  length  scale  of  the  probe  can  easily  approach  the  scale  of  the 
transverse  curvature.  As  mentioned  earlier,  the  determination  of  the  wall  shear  stress  is  also 
experimentally  difficult.  The  methods  that  were  used  are  noted  in  the  table. 

Three  methods  of  scaling  the  wall-normal  coordinate  have  been  commonly  used  for  the  mean 
velocity  profile  data.  The  first  method  is  to  use  the  planar  wail-normal  coordinate,  y+.  When  this  scaling 
is  used  for  small  5/a  the  results  are  very  similar  to  those  for  a  fiat  plate  boundary  layer/*2,41  although 
there  is  some  indication  that  the  coeflicients  of  *ht  planar  log  law  must  be  adjusted  to  lit  the  data.30,37 
As  the  curvature  ratio,  5/a,  increases,  the  logarithmic  relationship  for  the  mean  velocity  profile  remains 
irstaci.but  the  slope  of  the  logarilhmic  egion  deviates  from  the  planar  log  law3,1 3,12  as  shown  in  figure 
3.  To  quantify  this  effect,  Lueptow  et  al.10  compiled  dat"  from  several  sources  and  assumed  a  log  law 
of  the  same  'orm  as  the  planar  case,  i.e., 

Uj.0  (i/m)lny+  +  n.  (2-4) 

Thev  found  that  as  5/a— A ,  m  ap/  roaches  its  planar  value  of  K=0.4,  the  Von  Karman  coefficient. 
Likewise,  n  approaches  its  pianar  value  of  5.0  as  5/a— *1.  Both  coefficients  increase  with  increasing  5/a. 
In  tact,  in  is  a  linea:  lunctton  ot  6/a.  Tnis  linear  reiationshir  appears  io  break  down  ior  o/d>4G,  peihaps 
because  of  experimental  problems  when  measuring  praties  on  extremely  small  cylinders  or  because  ol 
reiaminarization. 

This  result  is  not  surprising  in  light  of  the  increased  coefficient  of  friction  measured  (and 
predicted,  as  seen  in  later  sections)  for  cylindrical  boundary  layers.  Recall  that  the  coefficient  ot  friction 
can  be  expressed  in  terms  of  the  friction  velocity  as  C(=*2  (U^U^2.  If  the  coefficient  of  friction  is  larger  in 
a  cylindrical  boundary  layer  than  in  a  planar  boundary  layer,  then  the  nondirnensional  free  stream 
velocity,  U^/U^,  must  be  smaller.  Unless  the  nondirnensional  cylindrical  boundary  layer  thickness,  S+,  is 
much  smaller  than  in  the  planar  case,  the  slope  of  the  tog  region  of  a  cylindrical  velocity  profile  must  be 
less  than  that  of  a  flat  plate  profile . 
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Table  1 ,  Measurements  of  the  Mean  Velocity  Profile  in  the  Cylindrical  Boundary  Layer 
(Values  for  8/a  and  x/a  were  estimated  based  on  information  presented  in  each  paper.  The  studies  listed 
are  in  order  of  decreasing  curvature  ratio,  8/a.  The  ratio  x/a,  where  x  is  the  axial  position  of  the 
measurement,  typically  decreases  as  the  curvature  ratio  increases  because  of  experimental  constraints.) 
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Table  1.  (Cont'd)  Measurements  of  the  Mean  Velocity  Profile  in  the  Cylindrical  Boundary '  ayer 
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b)  Poor  velocity  probe  resolution  calls  the  validity  of  the  results  into  question,  especially  near  the  wall. 

c)  Lueptow  et  al.10  include  characterization  of  the  mean  velocity  profile  for  over  40  profiles. 

d)  Willmarth  et  al.12  report  includes  a  total  of  14  profiles. 

e)  Joseph  et  al.35  include  characterization  of  the  mean  velocity  profile  for  over  60  profiles. 
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The  mean  velocity  profile  data  described  above  allows  the  evaluation  of  the  scaling  parameters 
described  earlier.  Lueptow  et  al.13  vaned  5/a  by  making  measurements  at  several  locations  along  the 
length  of  a  cylinder  at  constant  Ra.  Since  the  slope  of  the  log  region  of  the  boundary  layer  profiles 
varied  even  with  constant  Ra,  Ra  seems  to  be  an  inappropriate  scaling  parameter  for  cylindrical 
boundary  layers. 

Rao34,42  suggested  that  the  wall-normal  coordinate  of  the  planar  boundary  layer  log  law  be 
replaced  by  the  corresponding  wall-normal  coordinate  of  the  sublayer  law  for  an  axisymmetric  boundary 
layer  given  in  equation  (2-3)  Thus  y+  in  the  planar  law  should  be  replaced  by  Y+=[a+ ln(r/a)j  leaving  all 
coefficients  of  the  planar  iog  law  intact  (K=0.4,  B=5).  This  results  in  the  velocity  profile  in  the  log  region 
of  a  cylindrical  boundary  layer  of 

U+  =*  (1/K)  ln[a+  ln(r/a)]  +  B.  (2-5) 

Considerable  controversy  followed  the  introduction  of  Rao’s  wall-normal  coordinate.  Y+. 
Undoubtedly,  the  introduction  of  this  ucale  reflects  the  geometry  of  the  cylindrical  boundary  layer.  Since 
Y+  -  a+ln[l  +(y/a)]  ~  a+(y/a)  -  y+  for  large  a,  the  flat  plate  logarithmic  velocity  profile  is  recovered  from 
(2-5).  Chase43-44  supported  the  Rao  hypothesis,  since  it  suggests  a  physically  meaningful  mixing 
length  that  increases  linearly  with  distance  from  the  wail  only  for  yss  ,  while  it  increases  logarithmically 
further  out  in  the  boundary  laye  .  On  the  other  hand,  Bradshaw  and  Patel45  point  out  that  "there  is  no 
reason  why  the  fully  turbulent  flow  should  respond  to ...  [transverse]  curvature  in  the  same  way  as  the 
sublayer ..."  For  this  reason  they  dismissed  the  success  of  Rao's  log  law  (2-5)  in  matching  experimental 
data  as  fortuitous. 

Rao  and  Keshavan1 1  measured  mean  velocity  profiles  corresponding  to  (2-5)  for  cylinder 
diameters  ranging  from  1/16  inch  (0A6  cm)  to  5.5  inches  (14  cm).  Unfortunately,  poo'  probe  re  olution 
calls  their  results  tor  small  cylinders  into  question.  This  problem  aside,  they  found  that  the  coefficients  of 
(2-5)  to  be  dependent  upon  a+  and  Ra.  Femholz  and  Podtschaske46  also  suggested  that  the 
coefficients  in  (2-5)  are  not  constant  based  on  an  analysis  of  the  data  obtained  by  Willmarth  et  al.12 
Instead  they  proposed  that  the  coefficients  are  functions  of  a  curvature  ratio,  0/a,  where  0  is  the 
momentum  thickness.  At  curvature  ratios  of  S/a=0(l),  Adomaitis36  and  Furoya  et  al.HU  found  that  the 
coefficients  of  (2-5)  are  the  planar  values,  K  =0.4  and  B=5. 
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A  third  scaling  for  the  wall-normal  coordinate,  {y+  [1  +(y/2a)]j,  was  proposed  by  Richmond7 
based  on  cylindrical  geometry  and  a  hypothesis  that  streamlines  are  lines  of  constant  mean  velocity. 
Substituting  this  scaling  for  y+  in  the  planar  log  law,  Richmond,7  Yasuhara  38  Willmarth  and  Yang,14 
and  Oisso'inette  and  Mellor39  found  that  data  for  the  cylindrical  mean  velocity  profile  matched  the 
planar  velocity  profile  quite  well  for  8/a<2.  However,  Richmond's  data  for  large  5/a  fell  well  below  the 
planar  profile  in  the  log  region.  Rao34  cast  doubt  on  the  validity  of  the  Richmond  wall-normal  scaling. 

As  a  result  Richmond's  wall-normal  scaling  has  been  discarded  in  favor  o*  the  traditional  planar 
wall-normal  scaling  or  Rao's  wall-normal  scaling. 

Other  measurements  of  the  mean  velocity  in  a  cylindrical  boundary  layer  should  be  mentioned 
for  completeness.  Kwon  and  Prevorsek25  measured  the  velocity  profile  surrounding  a  continuous  fiber 
moving  through  quiescent  air.  The  small  diameter  of  the  fiber  resulted  in  a  boundary  layer  that  was  at 
least  two  orders  of  magnitude  thicker  than  the  fiber  diameter,  but  only  measurement  in  the  outermost 
portion  of  tho  boundary  layer  could  be  made  because  of  the  size  of  the  velocity  probe.  Joseph  et  al.35 
measured  the  mean  velocity  profile  for  cylinders  in  both  water  and  wind  tunnels.  In  determining  the 
exponent  of  a  power  law  approximation  tc  the  mean  velocity  profile,  they  found  that  the  exponent  is 
dependent  upon  the  radius  of  the  cylinder  but,  curiously,  is  independent  of  the  fluid  viscosity.  Finally, 
Willmarth  et  al 33  measured  isovelocity  contouis  around  cylinders  at  small  angles  of  yaw.  The  boundary' 
layer  on  the  leeward  side  was  nearly  10  times  as  thick  as  the  boundary  layer  on  the  windward  side  tor  a 
yaw  angle  of  only  1°. 

In  summary,  it  appears  that  the  planar,  Rao,  and  Richmond  wall-normal  coordinates  all  result  in  a 
logarithmic  relationship  for  the  mean  velocity  profile  in  the  inner  region  of  a  cylindrical  boundary  layer, 
although  the  planar  and  Rao  coordinates  are  preferred,  however,  no  wall-normal  scaling  collapses  data 
on  to  a  single  curve  for  all  transverse  curvatures.  Untortunately,  much  of  the  data  suffers  poor 
measurement  techniques  or  small  curvature  ratios  so  that  it  is  difficult  to  make  conclusions  on  the 
character  or  scaling  of  tire  mean  velocity  profile  for  cylindrical  boundary  layers. 

2.5  ANALYTIC  MODELS  OF  THE  VELOCITY  PROFILE  IN  THE  INNER  REGION 

Analytic  approaches  to  the  mean  velocity  profile  in  the  inner  region  of  a  cylindrical  boundary 
layer  often  depend  upon  the  assumption  of  similarity  with  planar  boundary  layer.;,  so  the  validity  of  the 
method  depends  upon  the  suitability  of  using  planar  boundary  layer  characteristics  to  describe  a 
cylindrical  boundary  lay^r.  Most  approaches  have  been  based  upon  eddy  viscosity  or  mixing  length 
closure  schemes  to  account  for  the  Reynolds  stress  term  in  (2-2).  The  eddy  viscosity,  F,  is  defined  so 
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that  the  represe  tation  of  the  turbulent  shear  stress  is  similar  to  the  representation  of  the  viscous  shear 
stress.  The  mixing  length  approach  relates  ihe  turbulent  shear  stress  through  a  turbulent  length  scale. 
Im.  Using  these  closure  schemes  the  Reynolds  stress  is  expressed  as 

_  2 

-uv  -  e  3U/dy  -  lm  l3u/dy|  3u/3y.  (2-6) 

Of  course,  the  eddy  viscosity  and  mixing  length  are  related  by  e  =  lm2  |dU/dy|.  .dost  closure  schemes 
for  cylindrical  boundary  layers  have  been  presented  in  terms  of  a  mixing  length,  so  that  is  the  format 
used  to  compare  the  schemes  in  table  2. 


The  earliest  attempts  at  using  a  mixing  length  for  closure  were  based  on  the  direct  application 
of  the  mixing  length  concepts  used  tor  planar  boundary  layers.  Sparrov/  et  a!.47  used  an  expression  of 
mixing  length  based  on  Von  Karman's  similarity  hypothesis  using  a  local  characteristic  length  scale 
based  on  consecutive  wall-normal  derivatives  of  the  mean  velocity  profile. 


Other  early  attempts  used  a  simple  planar  mixing  length,  l+  -  1,^^/v  -  Ky+>  where  K  is  the 
Von  Kaman  constant.  In  this  representation  the  mixing  length  and,  hence,  the  turbulent  length  scales 
are  proportional  io  ine  uiSia iu6  mumi  wan.  \jiiioV5f\ii  aim  .jUtGukm4^  expanded  (2-f)  in  a  Maclaunn 
series  and  then  applied  this  mixing  length  to  achieve  a  complex  expression  for  the  mean  velocity  profile 
in  terms  of  the  boundary  layer  thickness  and  8/a.  Their  expression  was  general  enough  to  account  tor 
either  concave  or  convex  transverse  curvature.  Reid  and  Wilson3  also  used  this  simple  planar  mixing 
length  to  derive  an  expression  for  the  boundary  layer  velocity  profile  that  was  a  function  of  the  sublayer 
thickness  because  of  a  matching  condition  at  the  overlap  of  the  sublayer  and  the  log  region.  They 
extended  their  analysis  to  include  the  effect  of  surface  roughness  on  the  velocity  profile.  Bradshaw  and 
Patel45  carried  out  an  analysis  of  the  log  region  using  a  planar  mixing  length  along  with  equation  (2-1). 
The  r  resultant  logarithmic  velocity  protile  required  an  additive  coefficient  dependent  upon  a+- 
Matsui27  also  used  the  planar  mixing  length  to  model  the  turbulent  shear,  although  th  >  constant  k  was 
determined  empirically. 


The  simple  planar  mixing  length  has  been  modified  to  include  a  factor  accounting  for  the 
velocity  profile  in  the  sublayer.  This  allows  the  use  ot  the  logarithmic  velocity  profile  in  the  sublayer  as 
well  as  in  the  log  region  of  a  boundary  layer  velocity  profile.  Cebeci49  used  the  Van  Driest  sublayer 
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correction55  for  the  mixing  length,  while  Patel55  used  a  factor  suggested  by  Landweber  and  Poreh57 
to  modify  the  logarithmic  velocity  profile  in  the  sublayer. 


Table  2.  Proposals  for  the  Mixing  Length  for  a  Cylindrical  Boundary  Layer 


Mixing  Length. Expression _ 

Planar  Mixing  Length: 

l+  -  k  (dU47dy+)  /  (d2U+/dy+2)  k  =■  0.4  (planar  value) 

1+  ■  Ky+  K  -  0.4  (plans'-  value) 


*+  *  Ky+ 

l+  -  Ky+[1-  6"  (V^+)] 
l+  -  Ky+[tanh(y+2/  X2 )]°  5 


K  *  0.22  (empirical) 

K  *  0.4  (planar  value) 
X+  -  26  (planar  value) 
K  -  0.4  (planar  value) 

X+  m  63N3 


Mixing  Length  Based  on  Rao  Coordinate: 


l+  **  k  [a+in(r/a)]  (r/a)°-5 


k  »  0.4  (planar  value) 


l+  -  K  [a+ln(rIa)j  [l  -  O'  i  la+ln(r/a)]  ^+1]  k  -  0  4  (planar  value) 

»  26  (planar  value) 

l+  =  K  [ajn(r/a)]  [tanh(y^2/  X  2 )]°  5  k  =  0.4  (planar  value) 

X  *  63/^3 

T 


Other  Cylindrical  Mixing  Lengths: 
u  «  Ky+  [f-  e-  (y+A-+)]  (a/r)05 

L  -  (KyJ  tr+/(r+  +  Kcy+)) 

!+  -  [c  ?+r(d.J+/dY+)]°  3 


K  »  0.4  (planar  value) 

\  •  \(aJ 

k  •  0.4  (planar  value) 
c  «  1 .5  (empirical) 
c  0.0274  (empirical) 


l+  =  (Ky+  (1  +  (y+/2a+)]}  (V'5 


k  =  0.4  (planar  value) 


Reference 


Sparrow,  Eckert  &  Minkowycz47 

Ginevdkii  &  Solodkin48 
Reid  &  Wilson2 
Bradshaw  &  Patel45 

Matsui27 
Cebeci 49 

Patel 50 


White,  Lessman  &  Christoph51 

Rao42 

Cebeci52 


Denli  &  Landweber2 


Granville52 

Eickhoff54 

Lueptow,  Leehey&  Stellinger12 
Lueptow  &  Leehey55 
Granville52  based  on 
Richmond7 
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Of  course,  the  application  of  planar  mixing  lengths  to  a  cylindrical  ooundary  layer  are 
immediately  subject  to  question.  The  planar  mixing  length  assumption  requires  that  momentum  transfer 
occur  over  a  length  scale  that  is  proportional  to  the  distance  from  the  wall.  However,  as  suggested  by 
Luxton  et  al.8  and  confirmed  using  flow  visualization  by  Lueptow  and  Haritonidis,31  the  scale  of  flow 
structures  in  a  cylindrical  boundary  layer  can  be  much  greater  than  the  diameter  of  the  cylinder  for  small 
cylinders.  In  other  words,  the  wall  is  less  likely  to  control  the  length  scales  in  a  cylindrical  boundary  layer 
than  in  a  planar  boundary  layer.  As  a  result,  using  a  mixing  length  that  is  proportional  to  the  distance  from 
the  wall  is  probably  inappropriate  unless  S/a  is  small. 

Rao42  used  the  similarity  law  (2-5)  to  derive  a  mixing  length  shown  in  table  2.  White, 
Lessrnan,  and  Christoph51  note  that  this  fonn  of  the  mixing  length,  which  is  smaller  than  its  planar 
equivalent,  is  logical  since  "a  cylinder  has  less  ability  to  create  turbulent  shear  than  a  plane  surface." 
However,  as  discussed  earlier  with  regard  to  equation  (2-5),  this  approach  is  based  on  a  weak 
assumption  that  the  wall-normal  scaling  in  the  viscous  sublayer  should  carry  over  to  the  log  region. 


A  different  approach  results  from  the  direct  substitution  of  the  Rao  coordinate,  Y+  -[a+ 
ln(r/a)],  for  the  planar  wall-normal  coordinate,  y+1  in  the  mixing  length  scaling.  Cebeci52  repeated  his 
analysis  incorporating  this  substitution  in  a  mixing  length  with  a  Van  Driest  sublayer  modification  as 
suggested  by  White  (discussion  following  reference  49).  Denli  and  Landweber8  used  a  mixing  length 
similar  to  that  used  by  Patel50  except  that  the  Rao  coordinate  was  used  instead  of  y+.  However,  using 
the  Rao  coordinate  in  place  of  y+  in  the  mixing  length  expression  is  subject  to  question  for  two  reasons. 
First,  as  noted  earlier,  the  Rao  wall-normal  coordinate  derives  from  the  scaling  in  the  viscous  sublayer.  Its 
application  to  the  log  region  does  not  necessarily  follow.  Second,  like  the  planar  case,  the  Rao 
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has  less  control  over  the  length  scales  than  in  a  planar  boundary  layer. 


Two  other  cylindrical  mixing  length  expressions,  noted  in  table  2,  have  been  proposed. 
Starting  with  equation  (2  1)  and  assuming  a  logarithmic  velocity  profile,  Granville53  developed  an 
expression  for  a  mixing  length  that  accounts  for  the  transverse  curvature  of  the  cylinder.  This 
expression  was  coupled  with  the  Van  Driest  sublayer  modification  noting  that  the  sublayer  modification 
correction  coefficient,  X,  should  be  a  function  o!  a+  instead  of  constant.  However  the  validity  of  this 
mixing  length  may  be  in  doubt,  since  its  derivation  depends  on  the  assumption  of  the  existence  of  a 
logarithmic  relationship  between  U+  and  y+.  Eickhoff54  proposed  with  little  justification  or  rigor  that  a 
>..'.0  of  radii  including  an  empirical  constant  should  be  multiplied  times  the  usual  planar  mixing  length  to 


18 


TR  8389 


include  the  effec  of  transverse  curvature.  This  mixing  length  is  based  on  the  weak  assumption  that  as 
an  eddy  moves  away  from  the  cylinder's  wall,  a  lower  velocity  fluctuation  results  in  a  cylindrical  boundary 
layer  than  in  a  planar  boundary  layer. 

Mixing  length  expressions  derived  from  the  similarity  laws  proposed  by  Lueptow  et  al.13 
and  Richmond7  are  included  to  complete  table  2. 

•  olMILARITY  LAWS 


Similarity  laws  for  the  mean  veloc-ity  profile  in  the  inner  region  are  presented  in  table  3.  As 
noted  earlier,  the  four  similarity  laws  based  on  a  planar  mixing  length  probably  do  not  appropriately 
account  for  the  effect  of  transverse  curvature.  However,  comparison  of  these  similarity  laws  to 
experimental  data  suggests  that  they  show  the  expected  trends  for  transverse  curvature.1 1 i53 

Three  similarity  laws  included  in  table  3  are  based  on  Rao's  coordinate,  Y+=[a+  ln(r/a)].  Denli 
and  Landweber8  used  this  coordinate  in  a  mixing  length  to  derive  a  similarity  law.  Rao's  similarity  law34  is 
based  on  the  substitution  of  the  sublayer  coordinate  into  the  planar  log  law  (equation  2-5).  Chase43 
offered  a  slight  modification  of  this  similarity  law  by  adopting  a  planar  term  to  allow  for  smooth  transition 
between  Rao's  sublayer  law  (equation  2-3)  and  the  similarity  law.  For  reasons  outlined  earlier,  the  laws 
based  on  a  mixing  length  incorporating  Rao's  wall-normal  coordinate  are  subject  to  scrutiny. 

The  last  four  similarity  lav/s  included  in  table  3  are  not  based  on  mixing  length  derivations  or 
Ran  s  wall-normal  coordinate.  Lueptow  et  al.1 0  proposed  a  similarity  law  based  on  an  assumption  of 
'c;  .start!  eddy  viscosity  in  the  boundary  layer.  This  results  in  a  log  law  that  is  identical  to  the  planar  log 
law,  equation  (2-4).  As  mentioned  eaniur,  m  aiiu  ii  are  oepecdent  upon  the  ratio  u/a  and  are  empincally 
determined.  The  validity  ot  the  constant  eddy  viscosity  approach  as  a  first-order  approximation  was 
confirmed  by  direct  measurement  of  the  eddy  viscosity.1 88  Using  the  method  o(  matched  asymptotic 
expansions,  Afzal  and  Narasimha88,59  derived  an  expression  for  the  mean  velocity  profile  that  is  similar 
to  the  planar  log  law.  As  a+  becomes  small  and  5/a  become:  large,  the  coefficients  of  the  log  law  are 
dependent  upon  these  parameters.  Thus,  the  similarity  law  is  very  similar  to  that  proposed  by  Lueptov/ 
et  al.10  although  the  laws  were  derived  using  two  different  approaches. 
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Planar  Mixing  Length: 


U+=  (1/K)  n(y+)+8 


K  »  0.4;  B  =5 


Planar  law 


U+=  (1/k)  ln[Q(y/a)/  Q(ys/a)]+  a+ln[(a+ys)/a]  K  =  0.4 


Reid  &  Wilson3 


U+~  (1/K)  in{4y+/[l  +(r/a)°-5]2}+B 


U+=  (1/k)  ln[4a+Q(y/a)]+  B 


k  =»  0.4;  B  function  of  a+  Bradshaw  &  Patel45 


k  =  0.4;  B  function  of  a+  Patel50 


Mixing  Length  Based  on  Rao  Coordinate: 

U+=*  (1/k)  ln[a+(r/a)°-5ln(r/a)]+B 


Based  on  Rao  Coordinate: 

U+=  A  ln[a+ln(r/a)]+B 


K  “  0.4;  B  function  of  a+  Denli  &  Landweber3 


A,B  functions  of  Ra>  a+  Rao;34  Rao  &  Keshavan1 1 


U+=  (1/K)  ln[a+ln(r/a)-J]-i  B 


K  *  0.4;  J  *  B  <*  5 


Chase43 


Other: 


U+=  A  ln(y+)+8 
U+-  A  ln(y+)+B 
U*=  A  ln{y41+(y/2a)]}+B 


A,B  functions  of  8/a  Lueptow,  Leehey,  &  Stellinger1 0 


A,B  functionsof  a+ 
not  specified 


Afzal  &  Narasimha53’50 
Richmond7 


U+=  (1/k)  ln(y+)+B 


K  =  0.4;  B  function  of  Ra  Yu30 


The  similarity  laws  due  to  Richmond  and  Yu  are  included  to  complete  table  3.  As  discussed 
earlier,  Richmond’s  law7  is  unlikely  to  be  useful.  The  proposal  by  Yu30  is  based  on  dimensional 
arguments  and  a  narrow  range  of  data.  Over  a  broader  range  of  transverse  curvature  ratios  it  fails. 

2.7  MEAN  VELOCITY  PROFILE  IN  THE  OUTER  REGION 

Like  a  planar  boundary  layer,  the  mean  velocity  profile  of  a  cylindrical  boundary  layer  deviates 
from  a  logarithmic  profile  in  its  outer  region.  The  mean  velocity  in  the  ouier  region  is  usually  displayed  in 

velocity-defect  coordinates,  <U00-U)/U1;  versus  y/6.  Willmarth  et  al.12  presented  velocity-defect  profiles 

over  a  wide  range  of  transverse  curvatures  as  shown  in  figure  2.  Their  data  clearly  indicate  that  the 
velocity  profiles  do  not  collapse  onto  a  single  curve  as  they  do  in  the  planar  case.  As  the  cylinder 

becomes  smaller  the  velocity-defect  profile  becomes  fuller.  As  8/a  approaches  one,  the  velocity-defect 
profiles  approach  the  planar  velocity  defect  profiles.  From  dimensional  arguments  Willmauh  et  al. 

propose  that  the  velocity  detect  should  be  a  function  of  both  y/6  and  8/a. 

Several  outer  laws  have,  been  proposed  tor  a  cylindncal  boundary  layer  as  ennumeratea  in  taoie 
4.  Afzal  and  Narasimha58-59  used  the  method  of  matched  asymptotic  expansions  to  obtain  a  defect 
law  that  is  quite  similar  to  the  planar  outer  law  except  that  the  multiplicative  coefficient  is  dependent 

upon  a+  and  the  additive  coefficient  is  a  function  of  both  5/a  and  a+.  Of  course,  the  coefficients 

approach  the  planar  vaiues  as  the  effect  of  transverse  curvature  diminishes  [5/a  =-0(1)  and  a+  very 
large].  Yu33  developed  an  outer  law  based  on  a  length  scale  dependent  upon  the  friction  velocity  and 
the  kinematic  viscosity.  Although  experimental  data  were  shown  to  collapse  onto  a  single  curve,  the 
range  of  transverse  curvature  ratios  is  very  small,  so  no  conclusions  can  be  drawn.  Based  on  similarity  to 

i  n  axisymmetric  wake,  Rao  and  Keshevan11  proposed  plotting  the  velocity  defect  versus  r+=(a+y)Ux/v. 

Denli  and  Landweber8  used  a  mean-llow  equation  simplil.ed  by  using  the  free  stream  velocity  in  the 
convective  term  (Oseen's  approximation)  and  an  eddy  viscosity  model  that  is  a  function  of  the 
longitudinal  coordinate  for  the  shear  stress  in  the  outer  layer  to  derive  a  complex  outer  law  requiring  six 

empirical  constants.  The  resulting  velocity  profiles  match  the  data  of  Willmarth  et  al.12  lor  8/a  ranging 
from  2  to  16.  Adomaitis36  developed  an  outer  law  based  on  the  planar  outer  law  and  using  Rao’s 

coordinates  and  replacing  8  with  the  displacement  thickness,  5*.  This  formulation  appears  adequate  tor 
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a  narrow  range  cf  5/a.  Adomaitis  also  documented  the  effect  of  free  stream  turbi  lence  on  t.ie 
bound  ary  layer  profile. 

Like  a  planar  boundary  layer,  the  deviation  of  the  velocity  profile  in  the  outer  layer  is  thought  of 
in  terms  of modification  of  the  logarithmic  velocity  profile  known  as  a  wake  law.  Rao  and  Keshavan1 1 
identified  two  versions  of  a  wake  component  in  a  cylindrical  boundary  layer.  The  wake  is  termed 
"negative"  if  the  boundary  layer  profile  drops  be'ow  the  log  law  as  the  distance  from  the  wall  increases 
and  "positive '  if  the  boundary  layer  profile  rises  above  the  log  law.  The  positive  wake's  appearance  is 
similar  to  a  planar  wake.  Although  Rao  and  Keshavan  considered  a  log  law  of  the  form  of  (2-  5),  these 
wakes  are  observed  regardless  whether  the  data  are  plotted  as  a  function  of  the  planar  wall-normal 
cooroinate  or  Rao's  wall-normal  coordinate. 


For  minimal  transverse  curvature  (S/a<2),  the  wake  is  positive  and  the  appearance  of  the  mean 
velocity  profile  is  nearly  identical  to  the  planar  profile.32'37'39'41  Chin  et  al 37  proposed  a  wake  law 

that  is  very  similar  to  the  planar  wake  lew  (see  table  4).  For  slightly  larger  curvature  ratios  [S/a=0(2)  he 
positive  wake  is  net  evident.14-36  Rao  and  Keshavan1 1  found  that  as  the  boundary  layer  developed, 
the  wake  went  from  negative  to  positive  at  a  given  Ra.  They  went  on  to  propose  that  the  radia' 

coordinate  at  which  the  velocity  profile  deviates  from  the  logarithmic  profile  is  dependent  upon  Ra. 

Their  explanation  of  this  phenomenon  involves  the  relationship  between  the  outer  portion  of  an 
axisymmetric  boundary  layer  ano  properties  of  an  axisymmetric  wake. 

On  the  other  hand  Afzal  arid  Narasimha68  call  existence  of  these  wakes  into  question  by 
proposing  that  the  flows  where  they  exist  may  not  be  fully  developed.  The  data  of  Willmarth  et  al.12 
spans  a  wide  range  of  transverse  curvatures  and  shows  a  posiVvo  wake  tor  large  transverse  curvature 

(6/a<5)  and  a  slight  tendency  toward  a  negative  wake  as  the  transverse  curvature  decreases  ( 5/a>9). 
However  they  point  out  that  the  negative  wake  in  a  cylindrical  boundary  layer  is  not  as  obvious  as  the 
positive  wake  in  the  planar  case,  in  spite  of  the  expectation  that  a  cylindrical  boundary  layer  should  be 
more  wake-like  than  a  plar.ar  boundary  layer.  As  a  result,  the  existence  of  a  wake  region  in  a  cylindrical 
boundary  remains  unclear. 
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Table  4.  Proposals  for  Outer  Laws  for  a  Cylindrical  Boundary  Layer 


Outer  Law 

Coefficients 

Reference 

(  U„-  U)/Ux  -  (-1/K)  ln(y/8)+C 

K=0.4;  C  =2.F 

Planar  law 

(  U«-  U)/Ut  -  -A  ln(y/5)+C 

A  function  of  a  t.; 

C  function  of  8/a,a+ 

Afzal  &  Narasimha^®-^ 

(  UM-  U)/Ux  .  (-1/K)  ln[y+exp{-U0</KUx)]+C 

K=0.4;  C  function  of  Ra 

Y>j30 

{ U«-  U)/Ux  « -A  lr(r+)+C 

A,C  functions  of  Ra 

Rao  &  Keshavan1 1 

( IW  U)/  *  flRaix/aRa.  ,'/(5+a)] 

Deal:  &  Landweber0 

( Uoo*  li)/u^  °  (*  w\)  ii ma  ii \  uww  jr)}tC  K»C.4|  0  cf  5/s 

Adomaitis36 

( U)/Ux  -  -  A  ln(y/5)+  C  (1-  cos(rcy/5)] 

A  -1.95;  0-0.9 

Chin,  Hulsebos,  & 
Hunnicutt37 

2.8  RESULTS  OP  ANALYTIC  MODELS  OF  A  CYLINDRICAL  BOUNDARY  I  .AYER 

Several  models  of  a  cylindrical  boundary  layer  have  been  used  to  evaluate  the  effect  of 
transverse  curvature  on  characteristics  of  a  boundary  layer  such  as  the  coefficient  of  friction,  the 
momentum  thickness  arid  the  displacement  thickness.  The  procedure  begins  with  the  development  of 
a  mean  velocity  p.ofile  based  on  one  of  the  mixing  length  closure  schemes  described  earlier  or  based 
on  the  assumption  of  a  particular  form  for  the  mean  velocity  profile.  From  the  velocity  profile  the 
coefficient  ot  trie  on,  Cf,  can  be  evaluated  from  the  momentum  integral, 


23 


TR  8389 


The  model  mean  velocity  profile  can  also  be  used  to  evaluate  the  integral  thicknesses.  The 
definitions  for  the  integral  thicknesses  differ  from  those  for  a  planar  boundary  layer  because  of  the  effect  of 
transverse  curvature.60  The  displacement  thickness,  8*,  and  the  momentum  thickness,  8,  are  given  by 
equations  (2-8)  and  (2-9),  respectively. 


S*  [1  +  (0*/2a)]  =  Ja+S[  (l-jf  )  (ijdr,  (2-8) 

6|1  t(0/2a)]  ,  Ja+5  ('  -  tfcj  (-§)]*•  <2  9> 

Table  5  indicates  the  functional  dependence  of  the  coefficient  of  friction  and  the  integral 
thicknesses  predicted  using  various  models  of  the  cylindrical  boundary  layer.  The  analyses  are  grouped 
according  the  model  used  for  the  mean  velocity  profile.  The  table  also  indicates  if  representative  mean 
velocity  profiles  were  presented  and  if  they  were  compared  to  experimental  results. 

The  earliest  attempts  to  model  a  cylindrical  boundary  layer  used  a  power  law  formulation  for  the 
mean  velocity  profile.  Millikan,61*  Landweber,6^  Eckert, 2  and  Sakiadis4  used  a  t/7-power  law 
formulation  assuming  similarity  to  a  planar  boundary  layer.  Karha  70  and  Liu  and  Dai63  varied  the 
exponent  of  the  power  law  in  attempts  to  refine  the  mean  velocity  prof'lo  to  account  for  transverse 
curvature.  The  power  law  description  of  the  boundary  iayer  profile  predicts  that  transverse  curvature  will 
increase  the  skin  friction  compared  to  the  planar  case.  However,  b  ^ause  of  the  crudeness  of  the 
velocity  profile  models,  the  predictions  are  not  quantitatively  correct. 


TKa  rrvs^rtlr  hnpA/H  An  fhn  nlnnor  mivlnn  lonnrth  oro  ni  itnrnQtirallu  ci  ihinrt  to  crn  itinw 

I  I  ICS  QliaiJllU  lirjuuio  UUJW  vn  mv  P'W  "iiAiny  'vm^iii  V*I  v  uuivuiMitvwi;  vv/wjwvi  vvv»»>) 


because  of  the  premise  of  similarity  of  the  mixing  length  in  the  cylindrical  case  to  the  mixing  length  in  the 
planar  case.  In  spite  of  this  weakness,  the  analyses  of  Sparrow  et  al.47  and  Cebeci49  stand  out  for  their 


*  Millikan’s  model  was  the  first  of  several  that  have  been  based  on  momen’um  integral  methods. 
While  these  methods  have  generated  considerable  interest,69  a  necessary  assumption  tor  their 
application  is  that  5«a.  Since  this  report  is  generally  concerned  with  situations  where  8>a,  discussion  of 
these  methods  is  omitted. 
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estimates  of  the  skin  friction  over  a  wide  range  of  cylinder  diameters,  Ra,  and  axial  positions,  Rx.  These 
results  indicate  that  the  coefficient  of  friction  decreases  with  increasing  Ra  and  Rx.  Sparrow  et  al. 
explained  this  in  the  following  way:  "...  the  difference  in  the  flow  field  about  a  cylinder  and  that  about  a  flat 
plate  is  due  to  the  expansion  in  flow  area  encountered  with  increasing  distance  from  the  cylindrical  surface. 
When  the  boundary  layer  is  very  thin  relative  to  the  cylinder  radius,  then  there  is  a  correspondingly  small 
area  change  across  the  boundary  layer,  and  the  cylinder  behaves  like  a  flat  plate.  On  the  other  hand,  when 
the  boundary  layer  is  thick  compared  with  the  cylinder  radius,  there  is  a  substantial  area  change  and  a  larger 
difference  between  flat  plate  and  cylinder...  Therefore,  the  larger  the  length  Reynolds  number  Rex  [Rx], 
the  thicker  the  boundary  layer  and  the  greater  the  effect  of  cylindrical  geometry.  The  cylinder  Reynolds 
number  Re^  [Ra]  is  a  direct  measure  of  the  radius  of  the  cylinder  and  so,  the  smaller  the  Re^  [Ra],  the 
greater  the  curvature  effect." 

In  spite  of  use  of  a  planar  mixing  length  mode!  in  these  analyses,  the  results  show  characteristics 
that  correspond  to  experimental  results.  Representative  velocity  profiles  presented  by  Sparrow  et  al. 
shov  the  shallow  slope  of  the  log  region  of  the  mean  velocity  profile  as  the  cylinder  radius  becomes  small. 
The  log  regions  of  their  velocity  profiles  also  show  a  slight  downward  curvature  that  is  also  vaguely  evident 
in  the  experimental  data  of  wmmartn  et  ai.12  for  very  smaii  cyiinuers.  A  similar  trend  is  seen  in  the 
representative  velocity  profiles  presented  by  Patel.50 

Using  planar  mixing  length  closure  models,  Reid  and  Wilson3  and  Matsui27  lound  the  unlikely 
result  that  the  average  skin  friction  is  a  function  of  Ra  alone  with  no  dependence  upon  tho  length  of  the 
cylinder.  Ginovskii  and  Solodkin48  show  qualitatively  similar  results  over  a  narrower  range  of  parameters. 
Reid  and  Wilson  also  considered  the  effect  of  wall  roughness  on  tho  skin  friction. 

The  use  of  cylindrical  mixing  length  models  noted  in  table  5  would  be  expected  to  provide  a  more 
appropriate  model  compared  to  the  planar  mixing  length  models.  Although  several  of  the  authors  who 
used  this  type  of  analysis  compare  their  mean  velocity  profile  to  experimental  data,  only  Denli  and 
Landweber0  go  so  far  as  to  calculate  the  skin  Indion  and  the  intogra!  thicknesses.  Their  calculations  agree 
with  the  experimental  results  of  Willmarth  et  a!.,12  although  this  may  be  a  result  of  the  use  of  this 
experimental  data  to  determine  the  empirical  coefficients  used  in  the  analysis. 
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Table  5:  Functional  Dependence  of  the  Coefficient  of  Friction  and  the 
Integral  Thicknesses  for  a  Cylindrical  Boundary  Layer 


BHS3X7sa9Beii£Ease=su3gcsBa: 

Mean  Velocity 

£Q£iIi£je 

HU1  integral 

fieieiencfl 

Power  Law: 

Emlil£a 

Friction13 

Thickness13 

Notes 

Millikan61 

- 

Rv 

Rv,x 

Rv  based  on  volume 

Eckert2 

- 

5/a 

8/a 

compressible 

Liu  &  Dai63 

- 

X 

X 

- 

Sakiadis4 

- 

- 

% 

- 

Planar  Mixing  Length  Closure: 

Cebeci49 

comparison 

Ra>  Rx 

- 

heat  transfer 

Ginevskii  &  Solodkin4® 

sample  prof. 

x/a,  Ra 

8/a 

- 

Matsui2? 

- 

Ra 

- 

extruding 

Patel59 

comparison 

- 

- 

Reid  &  Wilson3 

sample  prof. 

Ra 

- 

roughness 

Sparrow  et  a! 47  sample  prof. 

Cylindrical  Mixing  Length  Closure: 

Ra'  Rx 

heat  transfer 

Cebeci52 

comparison 

- 

- 

- 

UUIlh  i»V4I<’41(Vmwi 

rorTtnaricnn 

Vw*  ■  '«»»■  • 

- 

r- 

- 

Eickhoff54 

comparison 

- 

- 

Granville53 

- 

- 

- 

kinetic  energy 

Rao34 

comparison 

- 

- 

- 

Assumed  Velocity  Profile: 

Ackroyd64 

- 

Ra>  Rx 

- 

extruding 

White65 

comparison 

Ra'  Rx 

x/a,  Rx 

* 

White  et  al.51 

- 

Ra>  Rx 

- 

compressible 

Other  Methods: 

Abdelhaiim  et  al.*35 

- 

- 

- 

- 

Shanebrook  &  Sumner6? 

comparison 

- 

X 

- 

Yu30 

comparison 

Ra,  8/a 

Ra. 

- 

Notes: 

a)  "Comparison"  means  that  the  mean  velocity  profile  was  compared  to  experimental  data.  "Sample  prof." 
means  that  sample  mean  velocity  profiles  wore  provided  but  not  compared  to  experiments. 

b)  Parameters  shown  indicate  the  functional  dependence  of  the  coelficier.t  of  friction  or  the  integral 
thicknesses  (displacement,  momentum,  or  boundary  layer)  that  is  predicted  by  the  analysis. 
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Other  approaches  based  on  assuming  a  velocity  profile  instead  of  using  a  mixing  length  have 
boon  i"iCv  for  analysis  of  cylindrical  boundary  layer.  White65  substituted  Rao's  wall-normal  coordinate 
in  place  *.*  y ,  in  Spalding's  expression71  for  the  velocity  profile  in  a  planar  boundary  layer .  Ackroyd64 
corrected  errors  in  White’s  analysis  and  followed  a  similar  approach  for  the  problem  of  an  extruding 
cylinder.  The  results  for  a  cylinder  in  a  uniform  stream  and  for  an  extruding  cylinder  are  very  similar.  Both 
White  and  Ackroyd  prasent  results  over  a  wide  range  of  parameters  showing  that  the  coefficient  of 
friction  decreases  with  increasing  Ra  and  Rx  (see  figure  4) .  Unfortunately,  the  analysis  is  dependent 
upon  the  presumed  similarity  of  Spalding's  mean  velocity  profile  using  the  Rao  coordinate  with  the 
velocity  profile  in  a  cylindrical  boundary  layer.  Although  Ackroyd  clearly  shows  that  the  assumed  profile 

is  quite  different  from  experimental  data  for  large  6/a,  his  results  are  shown  in  figure  4  since  they  show 
the  trends  also  found  by  several  other  researchers.47'4^  White  et  al 51  carried  the  analysis  further  to 
include  compressibility. 


Finally,  three  other  methods  of  modeling  a  cylindrical  boundary  layer  were  used.  An 
entrainment  theory  based  on  fhe  momentum  equation  and  an  entrainment  equation  derived  from 
continuity  was  adapted  by  Shanebrook  and  Sumner®7  for  application  to  axisymmetric  boundary  lavers. 
GranvillQ68exte,?Jev.<  Uvs  work  substantially  by  incorporating  more  sophisticated  shape  parameter 
functions  and  velocity  profiles.  Yu30  developed  a  log  law  and  an  outer  law  similar  to  those  for  a  flat  plate 
and  used  these  laws  to  predict  the  boundary  layer  thickness  development  and  the  skin  friction  for 
cylinders  with  largo  Ra.  ^bdelhalim  et  al.66  developed  a  method  for  analyzing  the  flow  past  semi-infinite 
axisymmeti  o  iiiid  planar  bodies  with  a  blunt  leading  edge  condition  using  conformal  coordinates  and  the 


k-e  turouhnee  rrtodd  [i  tv/o-yqustion  closure  scheme).  Although  results  were  presented  for  the  planar 
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Figure  2.  Representative  Velocity  -Defect  Profiles  Showing  the  Influence  of  Transverse  Curvature 

[Willmarth  et  al.12 - Flat  plate.1'3  (U.  is  the  friction  velocity.) 

Reprinted  with  the  permission  of  Cambridge  University  Press,  copyright  1976.] 
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Figure  3.  Representative  Velocity  Profiles  in  Wall  Coordinates  Showing  the 
Influence  of  Transverse  Curvature 

[Willmarth  et  al.12  -  Flat  plate;1® . Sublayer  mean  velocity  profile  (equation  2-3  with 

a+=33.4).  (U*  is  the  friction  velocity.)  Reprinted  with  the  permission  of 
Cambridge  University  Press,  copynght  1976  ] 
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3.  TURBULENCE  IN  A  CYLINDRICAL  BOUNDARY  LAYER 


Although  the  emphasis  in  research  on  cylindrical  boundary  layers  has  centered  on  the 
measurement  and  analysis  of  the  mean  velocity  profile,  some  experimental  work  has  been  done  on  the 
turbulent  character  of  a  cylindrical  boundary'  layer.  Most  of  the  research  has  centered  on  the 
measurement  of  turbulence  in  a  cylindrical  boundary  layer  and  the  comparison  of  these  measurements 
with  a  planar  boundary  layer. 

3.1  REYNOLDS  STRESS 

Surprisingly  few  measurements  of  the  Reynolds  stress,  uv,  have  been  made  in  a  cylindrical 
boundary  layer.  Measurements  of  the  velocity  fluctuations  are  difficult  for  two  reasons.  First,  the 
measurement  ot  the  wall-normal  fluctuations  requires  the  use  of  X-  or  V-shaped  hot  wire  probes  that  are 
difficult  to  calibrate  and  use.  Second,  since  the  wall-normal  gradients  are  large,  the  measured 
wall-normal  velocity  is  strongly  affected  by  the  averaging  of  the  velocity  over  the  wall-normal  dimension 
of  the  hot-wire. 

The  Reynolds  stress  profile  in  a  cylindrical  boundary  layer  is  very  similar  to  that  in  a  planar 
boundary  layer  when  the  effects  of  curvature  are  small,  &'a<2.32>41  -72  As  6/a  gets  large,  though,  the 
Reynolds  stress  profile  becomes  quite  different  from  that  of  a  Hat  plate  as  shown  in  figure  5.10  The 
Reynolds  stress  drops  off  much  more  quickly  with  distance  from  the  wall  in  a  cylindrical  boundary  layer 
than  in  a  planar  boundary  layer.  Lueptow  et  al.  attribute  this  to  the  cylindrical  geometry.10  From 
equation  (2-1)  and  (2-2)  they  suggest  that  the  Reynolds  stress  should  be  a  (unction  of  1/r.  Plots  of  the 
Reynolds  stress  nondimensionalized  with  the  tree-stream  velocity  versus  a/r  appear  to  result  in  weak 
collapse  of  the  data  onto  a  single  curve  that  may  bo  considered  to  be  a  straight  line  as  a  first 
approximation.  Assuming  inertial  effects  are  small  throughout  the  boundary  layer,  equation  (2-1) 
suggests  that  the  slope  of  this  line  is  equal  to  the  coefficient  of  friction.  Using  this  method 
underestimates  the  wall  shear  stress,  but  the  degree  of  agreement  of  the  estimates  ot  the  coefficient  of 
friction  with  that  of  other  methods  is  quite  surprising  considering  the  necessary  assumptions.31 

The  Reynolds  stress  nondimensionalized  with  urms  and  v^j,  where  urms  ■=  V(uu),  is  nearly 
constant  across  the  entire  boundary  layer,  except  near  the  wall  and  near  the  outer  edge  of  the  boundary 


TR  8389 


iayer  where  it  necessarily  drops  off,1  °-32  Like  a  planar  boundary  layer  the  value  of  uv/(urmsvrms)  =0.5 
throughout  most  of  the  boundary  layer, 

3.2  TURBULENT  KINETIC  ENERGY 

The  kinetic  energy  of  the  turbulent  velocity  fluctuations  is  given  by 

Efc  =  (1/2)  p  UjUj  =  (1/2)p(uu+  w  +  ww  ),  (3-1) 

where  u,  v,  and  w  are  the  velocity  fluctuations  in  the  streamwise,  wall-normal,  and  spanwise  directions, 
respectively.  L:ke  the  measurement  of  the  wall-normal  velocity  fluctuations,  the  measurement  of  the 
spanwise  flue  nations  is  difficult. 

The  velocity  fluctuations  in  a  cylindrical  bour-.  ary  layer  are  nearly  identical  to  the 
measurements  cf  these  quantities  in  a  planar  boundary  layer  for  5/a<1 .41  >?2  For  larger  values  of  &a  the 
streamwise  velocsiy  fluctuation,  which  are  the  largest  contributor  to  the  turbulent  kinetic  energy,  appear 
to  drop  off  more  quickly  with  distance  from  the  wall  than  the  planar  case.32  On  the  other  hand,  the 
general  character  of  the  streamwise  and  wall-normal  velocity  fluctuations  near  the  wall  is  similar.to  that  of  a 
flat  plate  as  shown  in  figure  6.31  The  data  indicates  Reynolds  number  similarity  near  the  wall  when  the 
streamwise  velocity  fluctuations  nondimensionalized  by  the  Iriction  velocity  are  plotted  as  a  function  ot 
y +.  The  streamwise  velocity  fluctuations  rise  quickly  from  zero  at  the  wall  reaching  a  maximum  value  of 

urms/ux  =3-2  a1  y+  "13  3,1(1  dropping  oft  gradually  with  distance  from  the  wall  from  that  point  to  the  outer 
edge  of  the  boundary  layer.  Luxton  et  al.9  and  Lueptow  et  al.10  measured  a  slightly  smaller  maximum 
value  of  u,-ms/Ux,  although  the  discrepancy  may  be  related  to  probe  resolution  and  measurement 

techniques.  The  wall-normal  fluctuations  reach  a  maximum  value  of  v^g/U^.  -1 .  The  maximum  velocity 

fluctuations  and  the  distance  from  the  wall  where  they  occur  are  similar  to  channel  flow  and  planar 
boundary  layer  flow.  The  similarity  between  the  distribution  of  the  streamwise  and  wall-normal  velocity 
fluctuations  in  a  cylindrical  boundary  layer  and  other  wall-bounded  flows  suggests  that  the  mechanism 
for  the  generation  of  the  turbulence  at  the  wall  may  be  the  same. 
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Adomaitis36  found  that  the  magnitude  of  streamwise  velocity  fluctuations  increase  even  very 
close  to  the  wall  when  the  free  stream  turbulence  intensity  increases.  No  measurements  of  the 

spanwise  velocity  fluctuations  have  been  made  for  S/a>l. 

3.3  HIGHER  ORDER  VELOCITY  STATISTICS 

The  streamwise  skewness .TP/u^g3,  is  a  measure  of  the  direction  of  excursions  from  the  mean 
(positive  excursions  imply  positive  skewness,  negative  excursions  imply  negative  skewness,  and  a 
Gaussian  distribution  has  zero  skewness).  The  skewness  in  a  cylindrical  boundary  layer  is  similar  to  that  of 
other  wall-bounded  flows  with  a  value  between  0  and  -1  through  most  of  the  boundary  layer.9,74  Near  the 
wall  the  streamwise  skewness  becomes  positive  suggesting  the  presence  of  sweep  structures 
subsequent  to  a  burst  (see  section  4).  In  the  outer  part  of  the  boundary  layer  the  skewness  becomes 
more  negative  because  of  the  intermittency  of  the  turbulence.  The  wall-normal  skewness  is  slightly 
positive  throughout  the  boundary  layer74  indicating  that  the  largest  excursions  from  the  mean  are 
positive.  This  rnay  be  related  to  the  liftup  of  a  burst  near  the  wall  and  the  mild  positive  wall-normal  flow  in 
the  intermittent  region  near  the  edge  of  the  boundary  layer. 

The  streamwise  flatness,  u4/urms4,  is  a  measure  of  the  magnitude  of  excursions  from  the  mean 
(a  Gaussian  distribution  has  a  flatness  of  3;  flatness  greater  than  3  indicates  a  probability  distribution  with 
long  tails;  flatness  less  than  3  indicates  a  probability  distribution  with  short  tails  ).  The  streamwise  and 
wall-normal  flatness  are  nearly  Gaussian  throughout  most  of  the  boundary  layer.  Near  the  wall  the  flatness 
is  greater  than  3,  probably  because  of  the  burst-sweep  cycle.  In  the  outer  portion  of  the  boundary  layer 
the  flatness  is  greater  than  3  as  a  result  of  intermittency.9'74 

3.4  VELOCITY  SPECTRA  AND  AUTOCORRELATIONS 

The  spectrum  of  the  streamwise  velocity  fluctuations  is  indicative  of  the  energy  content  of  the 
velocity  fluctuations  at  a  particular  frequency,  f.  For  small  5/a  the  spectra  appear  to  be  similar  in  character  to 
the  planar  streamwise  velocity  spectra,32,41  although  it  is  difficult  to  make  direct  comparisons  because  of 
different  scaling  techniques.  Some  data  suggest  that  the  energy  in  a  cylindrical  boundary  layer  is  shifted 
to  higher  frequencies  compared  to  a  planar  boundary  layer.32 

Lueptow  ang  Haritonidis31  presented  the  energy  density  function,  E,  of  the  velocity  fluctuations 
such  that  the  area  under  the  spectrum  for  a  given  frequency  range  is  equivalent  to  the  turbulent  energy  in 
that  frequency  range.  The  ordinate  is  normalized  so  that  the  total  area  under  the  curve  is  1 ,  i.e.. 


33 


TR  8389 


f0  (fE/U^g2)  d(lnf)-1.  (3-2) 


This  format  allows  the  easy  identification  of  the  frequency  range  that  has  the  largest  contribution  to  the 
turbulent  energy.  Using  this  format,  Lueptow  and  Haritonidis  found  that  the  maximum  energy  content  of 
the  velocity  fluctuations  occurs  in  the  same  frequency  band  regardless  of  distance  from  the  wall  for  8/a=7 
as  indicated  in  figure  7.  This  suggests  that  turbulent  eddies  are  about  the  same  size  regardless  of  the 
distance  from  the  wall  supporting  a  constant  eddy  viscosity  model  for  the  turbulent  shear  stresses. 

Lueptow  and  Haritonidis31  found  that  an  outer  scaling  for  the  frequency,  t0=  fS/U^,  results  in  a 
better  collapse  of  streamwise  velocity  spectra  for  different  Reynolds  numbers  than  an  inner  scaling,  f+= 
fv/Ux2,  especially  at  low  frequencies  as  shown  in  figure  8.  Very  near  the  wall,  y+<lO,  they  found  that  an 
inner  scaling  worked  best.  On  the  other  hand,  Luxton  et  a!.9  found  no  Reynolds  number  similarity  when 
the  spectra  are  scaled  using  the  outer  scales  of  the  momentum  thickness  and  U^.  As  shown  in  figure  8, 
the  energy  content  of  the  wall-normal  velocity  fluctuations  is  substantially  less  than  the  streamwise 
velocity  tluctuations  and  the  maximum  energy  density  occurs  at  a  higher  frequency.31  The  two  spectra 
merge  at  high  frequencies. 


The  streamwise  velocity  autocorrelation  indicates  that  the  turbulent  integral  length  scales  in  the 
outer  portion  of  a  cylindrical  boundary  layer  are  less  than  half  the  size  of  the  length  scales  in  a  planar 
boundary  layer.32  The  data  of  Luxton  et  al.9  suggest  that  length  scales  derived  tram  the  autocorelation 
appear  to  be  slightly  dependent  upon  the  location  in  the  boundary  layer.  They  are  about  twice  as  large  at 
y+-540  than  they  are  at  y+=lO  for  8/a»27.  Comparing  these  two  studies  it  appears  that  the  length  scales, 
L,  are  much  larger  for  large  8/a.  At  y/5=0.5,  the  autocorrelation  falls  to  1/e  at  1=0.088  for  S/a=i  .6,  but 
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scales  is  minor  for  large  8/a. 


3.5  WALL  PRESSURE 


The  wall  pressure  spectrum,  autocorrelation  and  convection  velocity  for  boundary  layers  with 
slight  transverse  curvature  (8/a<0.1)  are  in  agreement  with  data  obtained  for  planar  boundary  layers  and 
pipe  flow75.76  For  larger  8/a ,  Willmarth  and  Yang14  and  Willmarlh  et  al.12  investigated  the  wall  pressure 
fluctuations  using  an  array  of  0.06-inch  (0.15  cm)  flush-mounted  piezoelectric  pressure  transducers  on  a 
3-inch  (7.G  cm)  diameter  cylinder  (S/a=2)  and  on  a  1-inch  (2.5  cm)  diameter  cylinder  (  S/a=4).  The 
convection  velocity  of  wall  pressure  fluctuations  derived  from  the  space-time  wall  pressure  correlation  was 


34 


TR  8389 


found  to  be  nearly  the  same  as  that  for  a  planar  boundary  layer  increasing  with  streamwise  separation  from 
O.eiJ^  to  0.8Uto.  However,  since  the  mean  velocity  profile  converting  the  eddies  is  fuller  because  of  the 
transverse  curvature,  the  turbulent  eddies  must  be  smaller  to  result  in  the  same  convection  velocity  as  in  a 
planar  boundary  layer.  Consistent  with  this  idea,  the  power  spectra  of  the  wall  press,  're  fluctuations  for  a 
cylindrical  boundary  layer  contain  a  greater  energy  density  at  higher  frequencies  than  a  planar  boundary 
layer  as  shown  in  figure  9. 

Constant  wall  pressure  correlation  contours  in  a  cylindrical  boundary  layer  are  nearly  circular  in  the 
streamwise-spanvvise  plane.14  This  is  in  contrast  to  the  elongation  of  the  contours  in  the  spanwise 
direction  measured  in  a  flat  plate  boundary  layer  as  shr  wn  in  figure  10.  Willmarth  and  Yang  suggest  the 
following  explanation.  Consider  a  large  eddy  adjacent  to  the  curved  wall  of  a  cylhder.  The  mean  velocity  at 
the  spartwise  sides  of  the  eddy  is  necessarily  larger  in  a  cylindrical  boundary  layer  than  in  a  planar  boundary 
layer,  since  the  distance  from  the  wall  to  the  side  of  the  eddy  is  greater.  This  results  in  a  shearing  motion 
on  the  sides  of  the  tddy  and  reduces  its  iransverse  scale.  However,  the  data  of  Willmarth  et  al.12 
confuses  the  issue  somewhat.  They  measured  wall  pressure  contours  for  8/a~4  that  were  more 
elongated  in  the  sparrwise  direction  than  measured  by  Willmarth  and  Yang  for  5/3=2  but  less  than  in  the 
planar  case  ( 8/a»0). 


The  difference  in  the  wall  pressure  correlation  contours  between  a  cylindrical  boundary  layer  and 
a  planar  boundary  layer  suggest  differences  in  the  wavenumber  spectra  of  lha  wall  pressure  for  the  two 
cases.  As  noted  by  Blake80  with  regard  to  the  expression  for  the  wavenumber  spectrum  of  the 
mean-shear-turbulence  interaction  term,  the  spectrum  of  pressure  fluctuations  in  the  streamwise  direction 
is  enhanced  relative  to  the  spanwise  direction  in  the  planar  case.  This  results  in  spatial  pressure 
correlations  that  reflect  smaller  scales  in  the  streamwise  direction  than  in  the  spanwise  direction.  The  larger 


scales  In  the  sparrwiss  direction  arc  evident  in  the  elongation  of  the  contours  in  the  spanwise  direction. 


Although,  an  extension  of  this  approach  to  a  cylindrical  boundary  layer  is  beyond  the  scope  of  this  work.,  it 


seems  likely  that  the  differences  in  the  wall  pressure  correlation  contours  displayed  in  figure  10  are  related 
to  differences  in  the  wavenumber  spectra  for  the  wall  pressure. 


Two  analyses  related  to  the  low-wavenumber  spectrum  of  the  wall  pressure  have  been 
undertaken.  Chase  and  Noiseux81  related  the  turbulent,  wall  pressure  fluctuations  in  a  cylindrical 
boundary  layer  to  nonlinear  fluctuating  velocity  products,  or  Reynolds  stress  sources.  Using  a 
perturbation  approach  they  derived  integral  relations  expressing  the  wall  pressure  amplitude  at  c  given 
wavenumber  and  frequency  as  an  integral  over  the  nonlinear  sources.  The  expressions  were  developed 
for  the  low-wavenumbur  domain  as  expansions  in  a  parameter  U^k/G)"  1 ,  where  k  is  the  streamwise 
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wavenumber  and  o)  is  tlte  angular  frequency.  Although  this  work  provides  a  analytic  framework,  its 
extension  requires  the  modeling  of  the  source  spectra.  Dhanak82  investigated  the  low-wavenumber 
domain  of  pressure  fluctuations  in  a  cylindrical  boundary  layer  for  8/a«1 .  Beginning  with  the  Lighthill 
formulation  in  cylindrical  coordinates,  the  low-wavenumber  approximation  to  spectral  density  was  found  in 
terms  of  a  product  uf  source  terms  and  a  response  function  related  to  the  cylindrical  geometry.  Unknown 
functions  of  frequency  appearing  in  the  analysis  will  need  to  be  determined  experimentally  to  advance  this 
work, 


Willmarth  et  al.33  investigated  the  wall  pressure  on  cylinder  subject  to  slight  cross-flow.  For  a  yaw 
an^.le  of  2,4°  they  found  that  the  root-mean-square  wall  pressure  is  15%  higher  on  the  windward  side  and 
35%  lower  on  the  leeward  siae  than  on  a  cylinder  with  zero  yaw.  The  windward  side  wall  pressure 
contained  more  energy  in  the  high  frequency  components  than  the  leeward  side.  Interestingly,  periodic 
oscillations  of  the  wali  pressure  were  not  detected  indicating  an  absence  of  vortex  shedding  for  small  yaw 
angles. 

3.6  FLUCTUATING  WALL  SHEAR  STRESS 

The  fluctuating  wall  shear  stress  has  not  been  measured  for  a  cylindrical  boundary  layer.  This  is 
largely  due  to  the  difficulty  in  the  application  of  the  usual  methods  for  the  measurement  of  wali  shear  stress 
to  small  cylinders. 

3.7  INTERM  RTENCY 

The  intermitfency  at  the  outer  edge  of  a  cylindrical  boundary  layer  was  measured  by  Lueptow  and 
Harrtonidis0 1  for  8/a~7.  Using  a  detection  scheme  based  on  the  square  of  the  first  derivative  of  the 
velocity,  they  found  that  the  instantaneous  location  or  the  interface  between  turbulent  and  non-turbulant 
flow  is  Gaussian  in  character  just  like  a  planar  boundary  layer  and  that  the  data  collapse  onto  a  single  curve 
fur  several  Reynolds  numbers.  Howevc:,  the  mean  location  of  the  interface  was  at  y/S-1 .0  compared  with 
y/8-0.8  for  a  flat  plate  boundary  layer.  They  ar.:  bute  this  difference  to  more  energetic  streamwise  velocity 
fluctuations  further  out  in  a  cylindrical  boundary  layer  than  in  a  planar  boundary  layer  and  to  the  "filling  out" 
of  'he  boundary  layer  with  turbulent  eddies  since  'he  cy.inder  does  not  constrain  the  motion  of  the  eddies 
i;  e  a  flat  plate  does. 
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3.8  RELAViiNARIZATiON 


Substantial  discussion  regarding  the  relaminarization  ot  a  turbulent  cylindrical  boundary  layer  has 
resulted  from  consideration  of  extremely  small  cylinders  where  5/a— and  Ra— >0.  Cebeci40  tried  to 
relate  th»  minimum  Reynolds  number  based  on  momentum  thickness  that  is  a  necessary,  though  not 
sufficient,  condition  for  turbulence  to  this  problem.  For  Ra=10,  he  predicted  that  the  displacement 
thickness  would  need  to  be  at  leas?  32  times  greater  than  the  cylinder  radius  for  the  flow  to  remain 
turbulent.  Based  on  extrapolation  of  their  experimental  results,  Rao  and  Keshavan1 1  proposed  that 
relaminarization  may  occur  when  Rg  *1 5,000  because  the  coefficient  of  friction  for  laminar  flow  will  match 
that  for  turbulent  flow  very  far  downstream.  This  critical  Reynolds  number  is  close  to  the  value  of 
Ra=l  1 ,000  below  which  smail  disturbance  will  be  damped  based  on  hydrodynamic  stability  of  a  cylindrical 
boundary  layer.83  White55  argued  that  similar  coefficients  of  friction  for  laminar  and  turbulent  flow 
suggest  relaminarization  at  Ra-3.  Turbulent  cylindrical  boundary  layers  have  been  measured  tor  Ra  as 
small  as  96  and  140.7,9 


Patel50  pointed  to  The  rather  remarkable  similarity  between  the  influence  of  a  favourable 
pressure  gradient  and  that  of  transverse  wall  curvature.’*  Using  a  shear  stress  gradient  parameter  reiated 
to  relaminarization  in  a  planar  houndary  layer  with  a  favorable  pressure  gradient,  he  predicted  that  a 
cylindrical  boundary  layer  would  relaminarize  for  a+<28.  Afzal  and  Narasimha58  supported  this  proposal 
based  upon  the  “inadequate  energy  supply  from  the  surface  [wall]."  However,  Willmarih  et  al.12  pointed 
out  that  the  apparent  similarity  in  the  mean  velocity  profiles  for  a  cylindrical  boundary  layer  ai  id  a  planar 
boundary  layer  with  a  favorable  pressure  gradient  result  from  quite  different  physical  phenomena. 
Transverse  curvature  causes  a  full  velocity  profile  in  a  cylindrical  boundary  layer  whereas  free  stream 
acceleration  causes  the  same  eflect  in  a  planar  boundary  layer  with  a  favorable  pressure  gradient.  In  any 
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Figure  5.  Profile  of  the  Reynolds  Stress  in  a  Cy.  ndrical  Boundary  Layer 

[Lueptow  et  al 1  ® - Flat  plate. ^  Reprinted  with  the  permission  of  the 

American  Institute  of  Physics,  copyright  1985.) 
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Figure  6.  Profile  of  the  Streamwise  (open  symbols)  and  Wall-Normal  Turbulence 
Intensity  (filled  symbols)  in  a  Cylindrical  Boundary  Layer 

[Lueptow  and  Haritonidis.31  A  R0=2OCiO;  Q  Rq~3300;  OR0=56OO. -  Flat  plat  71 

Reprinted  with  the  permission  of  the  American  Institute  of  Physics,  copyright  1987.] 
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Figure  7.  Streamwise  and  Wall-Normal  Velocity  Spectra  Near  the  Wall  in  a  Cylindrical 
Boundary  Layer:  -  Rg=2000,  y+»37; - Rg=3300, 


Boundary  Layer:  -  Rg=2Q00,  y+ss37; - Rg=3300, 

y+=32;  Rq-5600,  y+=47.  Wall-Normal  Velocity  Spectra:  Rg^2000,  y+~37 


[  Lueptow  and  Haritonidis.31  Reprinted  with  the  permission  of  the 
American  Institute  of  Physics,  copyright  1987.] 
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Figure  8.  Streamwice  Velocity  Spectra  at  Ditlerei;!  Distances  Irom  the  Wall  in  a  Cylindrical 

Boundary  Layer:  Rg«2000:  — — y+=2; - y^—13, 

y+-37; - y+~73, - y+~1 58 +  y+~399 

[Lueptow  and  Hartonidis.31  Reprinted  with  permission  of  the  American  Institute 

of  Physics,  copyright  1987.] 
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Figure  9.  Wall  Pressure  Spectrum  in  a  Cylindrical  Boundary  Layer 
[Wilimarth  and  Yang.1'1  Data  points  are  measured  for  a  cyiindricai  boundary  iayer  o/a~2. 

Solid  and  dashed  lines  are  two  different  measurements  for  a  planar  boundary  layer.77,78  <fi(co)  is  the 
power  spectrum,  co  is  the  angular  frequency,  and  p  is  the  turbulent  wall  pressure.  Reprinted  with  the 
permission  of  the  Cambridge  University  Press,  copyright  1970.] 


Figure  10.  Effect  of  Transverse  Curvature  on  Contours  ol  Constant  wall  Pressure  Correlation 
[W'llmarth  and  Yang.14  Solid  and  dashed  lines  are  two  different  measurements  for  a  planar  boundary 
layer.77,7^  Rpp(*i  .x3,0)  is  the  normalized  spatial  wall  pressure  correlation  at  zero  separation  time; 
x1  is  the  streamwise  coordinate;  X3  is  the  spanwise  coordinate.  Reprinted  with  the  permission  of 
the  Cambridge  University  Press,  copyright  1970.] 
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4.  S  TRUCTURE  OF  TURBULENCE  IN  A  CYLINDRICAL  BOUNDARY  LAYER 

Like  a  planar  boundary  layer,  an  understanding  of  the  underlying  structure  of  the  turbulence  in 
a  cylindrical  boundary  layer  is  fundamental,  yet  difficult.  Veiy  little  information  is  available  on  the 
structure  of  tudoulence  in  a  cylindrical  boundary  layer,  although  soma  aspects  are  becoming  more  clear. 

4.1  EVENT  DETECTION 

Lueptow  and  Haritonidis31  used  VITA,84  or  variable  interval  time  averaging,  an  event  sampling 
technique  based  on  the  detection  of  layers  of  high  shear,  to  detect  burst  cycles  which  are  believed  to 
be  responsible  for  generating  most  of  the  turbulent  velocity  fluctuations  in  a  planar  boundary  layer 85 
The  ensemble  averaged  events  detected  using  VITA  for  5/a=8  resemble  those  constructed  for  a  planar 
boundary  layer  and  channel  flow.  This  similarity  suggests  that  the  mechanism  for  the  generation  of 
tuibulence  at  the  wall  is  similar  in  all  wall-bounded  flows.  The  scaling  for  the  frequency  of  bursts 
detected  using  VITA  can  be  used  io  indicate  whether  the  mechanism  for  the  generation  of  turbulence 
scales  with  inner  or  outer  variables.  Altfwugh  their  data  are  inconclusive,  Lueptow  and  Haritonidis  found 
that  outer  scaling  provides  somewhat  better  collapse  than  inner  scaling,  suggesting  that  the  outer  flow 
affects  events  at  the  wall.  By  comparison,  the  burst  frequency  in  a  planar  boundary  layer  scales  with 
inner  variables. 

Lueptow  and  Haritonidis31  also  carried  out  uv-quadrant  event  detection  based  on  sorting  pairs 
of  u  and  v  into  appropriate  quadrants  of  the  u-v  plane  to  represent  contributions  to  the  Reynolds 
stress.86  Quadrants  2  and  4  are  associated  with  components  of  the  burst  cycle.  Like  VITA  detected 
events,  ensemble  averages  ol  uv-quadrant  detected  events  are  similar  to  those  of  other  wall-bounded 
fiows.  in  addition,  the  fractional  contribution  of  each  quadrant  in  the  u-v  plane  to  the  Reynolds  stress 
near  the  wall  (y+=39)  is  qualitatively  simile'  to  that  for  a  planar  boundary  layer.  Specifically,  quadrant  2 
events  (negative  u  and  positive  v)  represent  the  largest  contribution  to  the  Reynolds  stress  and 
correspond  to  the  lift-up  of  low  speed  fluid  during  the  burst  cycle.  Again  the  similarity  of  the  results  of 
event  detection  in  a  cylindrical  boundary  layer  using  the  jv-quadrant  technique  with  those  of  a  planar 
boundary  layer  suggest  that  the  mechanism  for  the  generation  of  turbulence  near  the  wall  is  similar. 

4.2  SPACE-TIME  VELOCITY  CORRELATIONS 

Measurements  of  the  correlation  of  velocity  at  different  azimuthal  locations  in  a  cylindrical 
boundary  layer  were  made  by  Lueptow  and  Haritonidis.31  In  these  experiments  the  streamwise  velocity 
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was  measured  using  three  hot-wire  probes  located  near  the  wall  at  three  a2imuthal  locations  that  were 
90°  apart  and  at  several  streamwise  spacings  for  flow  conditions  resulting  in  S/a=8.  The  results  of  the 
experiments  indicate  that  the  correlation  between  any  two  probes  is  negative.  This  result  is  curious 
since  the  spanwise  correlation  in  a  planar  boundary  layer  is  positive.  The  authors  suggest  that  the 
negative  correlation  may  be  a  result  of  the  splitting  of  a  sweep  of  tiuid  around  the  cylinder  pushing 
low-speed  fluid  ahead  of  ft  and  resulting  low-speed  and  high-speed  flow  90°  apart  from  each  other.  As 
would  be  expected,  velocities  at  probes  separated  by  90°  are  better  correlated  than  velocities  at  probes 
on  opposite  sides  of  the  cylinder  (separated  by  180°).  The  velocities  appear  correlated  over  distances 
of  at  least  13  diameters  suggesting  that  high-  and  low-speed  lumps  of  fluid  travel  in  adjacent  parallel 
paths.  The  advection  velocity  determined  from  the  cross-correlations  was  nearly  the  same  as  the  local 
mean  velocity. 

Finally,  VITA  event  detection  on  the  velocity  at  the  probes  positioned  near  the  wall  (y+=8)  and 
probes  positioned  further  from  the  wall  (y+*77)  was  used  to  determine  the  relationship  between  inner 
and  outer  events.  Regardless  of  whether  an  upstream  event  occurs  near  the  wall  or  a  short  distance 
from  the  wall,  it  seems  to  be  a  precursor  to  downstream  events  suggesting  an  interaction  between  the 
outer  flow  and  the  flow  near  the  wall. 

4.3  FLOW  VISUALIZATION 

Flow  visualization  of  a  cylindrical  boundary  layer  has  been  accomplished  by  moving  a  flexible 
cylinder  axially  through  a  quiescent  fluid.  Lueptow  and  Haritonidis31  used  a  system  in  which  a  very  long 
O-ring  moved  around  four  pulleys  located  at  the  comers  of  a  rectangular  frame.  One  of  the  pulleys  was 
driven  by  a  dc  motor.  The  frame  was  lowered  into  a  tank  of  water  so  that  the  bottom  two  pulleys  and  the 
portion  or  the  O-ring  uoiwooi'i  iheni  wtrrts  SubfiWQwu.  A  “hydrogen  bubble"  Wife  WS5  positioned 
perpendicular  to  and  just  above  the  O-ring  as  it  moved  through  the  tank  of  water  (  S/a~24).  In  spite  of 
the  buoyancy  of  the  bubbles,  groups  of  bubbles  were  occasionally  carried  by  the  turbulence  below  the 
O-ring  as  shown  in  figure  1 1 .  The  frequency  of  these  events  is  of  the  same  order  of  magnitude  as  the 
frequency  of  VITA  events  in  which  the  fluid  is  decelerating  in  the  streamwise  direction.  The  authors 
point  out  that  this  type  of  turbulent  transport  is  quite  different  from  that  of  a  planar  boundary  layer  where 
the  motion  of  large-scale  structures  is  constrained  by  the  wall.  As  8/a  becomes  large,  these  coherent 
structures  can  travel  from  one  side  of  the  cylinder  to  the  other  with  relative  ease. 

Other  researchers  have  used  birefringent  liquids  to  observe  the  flow  field  surrounding  moving 
filaments.  Gebart  and  Hornfeldt87  studied  annular  flow  in  a  pipe  with  a  nylon  filament  at  the  centerline  of 
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the  pipe  acting  as  a  moving  inner  boundary  condition.  The  relationship  of  their  results  to  a  cylindrical 
boundary  layer  are  difficult  to  ascertain  because  of  the  outer  pipe  wall  boundary  condition.  Sakiadis86 
provides  one  photo  of  the  flow  field  surrounding  a  moving  thread,  but  it  appears  that  the  flow  pictured  is 
laminar. 
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FlClirS  ^otrhoc  nf  tha  Mntirm  nf  I  arno  Rtmrturac  in  a  f!\/linHrir.al  RminHarv  I  avor 

[Upper:  Side  view  of  a  cylinder  (a  long  O-ring)  moving  to  the  right  through  quiescent  water.  A  wire  just 
above  the  cylinder  and  perpendicular  to  the  plane  of  the  paper  generates  bubbles.  Lower:  End  view  of 
a  cylinder  (O-ring)  moving  out  of  the  plane  of  the  paper.  A  wire  just  above  the  cylinder  and  in  the  plane 
of  the  paper  generates  bubbles.  In  both  cases  the  sketches  indicate  a  sequence  of  four  positions  of  a 
small  group  of  bubbles  as  the  turbulent  fluctuations  carry  them  below  the  cylinder  (positions  1  and  4  are 
labeled).  Note  that  gravity  is  downward  in  the  sketches,  so  the  turbulent  transport  carries  the  bubbles  in 
a  direction  opposite  to  their  buoyancy.  The  dashed  lines  indicates  the  approximate  edge  of  the 
boundary  layer.  Sketches  are  based  on  flow  visualization  by  Lueptow74  for  a+=30,  Rg*600,  8/a=20, 
time  between  positions=0,5  seconds.] 
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5.  SUMMARY 


A  description  of  the  boundary  layer  on  a  cylinder  in  axial  flow  has  evolved  over  the  past  30  years 
of  research  on  the  subject,  but  t:  ^  oresent  understanding  of  the  flow  lags  behind  that  of  other 
wall-bounded  flows.  This  is  a  consequence  of  the  difficulty  brought  about  by  the  introduction  of  an 
additional  length  scale,  the  transverse  curvature. 

In  the  viscous  sublayer,  the  transverse  curvature  results  in  an  equation  for  the  mean  velocity 
that  is  different  from  the  planar  law,  although  it  approaches  the  planar  law  as  the  radius  of  the  cylinder 
increases  (equation  2-3).  Away  from  the  wall  the  effect  of  transverse  curvature  becomes  pronounced  as 
the  mean  velocity  profile  becomes  fuller.  When  plotted  in  planar  wall  coordinates  this  appears  as  a  log 
law  relationship  with  the  slope  of  the  logarithmic  region  decreasing  as  the  transverse  curvature 
decreases.  Yet  the  parameter  of  interest,  the  transverse  curvature,  is  precisely  the  most  troublesome 
parameter  of  the  flow.  The  appropriate  scaling  for  the  wall-normal  coordinate  has  not  been  clearly 
defined. 


Like  scaling  laws,  appropriate  nondimensional  parameters  remain  elusive.  Although  several 
nondimensional  parameters  were  suggested  early  in  this  paper,  the  evidence  supporting  one  over  the 
others  is  weak.  The  radius  Reynolds  number ,  Ra,  has  already  been  shown  to  be  the  least  appropriate 
nondimensional  parameter,  since  it  omits  any  effect  due  to  varying  wall  shear  stresses  or  integral 
thicknesses  at  different  axial  positions.  For  this  reason  a+  appears  to  be  an  attractive  parameter,  but 
Lueptow  et  al.10  point  out  that  a  wide  range  of  mean  velocity  profile  data  collapse  better  with  8/a  than 
with  a+.  Both  8/a  and  a+  suffer  from  the  difficulty  that  measurements  of  the  boundary  must  be  made  to 
determine  these  parameters.  This  is  similar  to  the  problem  of  using  the  momentum  thickness  instead  of 
the  streamwise  position  for  the  length  scale  in  the  Reynolds  number  for  a  planar  boundary  layer  -  the 
momentum  thickness  length  scale  requires  making  a  measurement.  The  laminar  cylindrical  boundary 
layer  parameter,  overcomes  this  deficiency.  Except  for  Denli  and  landweber,0  employment  of  this 
parameter  has  been  minimal  so  its  usefulness  is  uncertain, 

More  certain  is  the  wake-like  character  of  a  boundary  layer  on  a  very  small  cylinder.  Denli  and 
Landweber8  describe  it  best.  "II  the  boundary-layer  thickness  relative  to  the  transverse  radius  of 
curvature  is  large,  the  cylinder  may  be  considered  as  a  small  vorticily-  and  turbulence-producing 
disturbance.  Consequently,  the  flow  might  be  considered  similar  to  a  wake  flow  with  the  modified  inner 
boundary  condition.  The  important  difference  is  that  the  drag  generating  the  wake  is  a  function  of  the 
longitudinal  coordinate..."  luxtcn  et  al.9  describe  a  cylindrical  boundary  layer  further  as  a  "continuously 
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regenerated  wake."  The  wake-like  character  o(  a  cylindrical  boundary  layer  is  clear,  but  the  existence  of 
an  cuter  wake  law  is  subject  to  debate  as  discussed  earlier. 

Although  the  outer  flow  is  wake-like,  the  wall  acts  to  continuously  convert  mean  flow  energy 
into  turbulent  energy  at  the  wall  of  the  cylinder.  In  fact,  a  cylindrical  boundary  layer  is  very  effective  at 
this  conversion  of  energy  as  evidenced  by  a  higher  coefficient  of  friction  for  a  cylindrical  boundary  layer 
than  for  a  planar  boundary  layer  at  the  same  Reynolds  number.  Detection  of  the  burst  cycle  (VITA  and 
uv-quadrant  detection  schemes)  as  well  as  measurements  of  the  turbulence  intensity  indicate  that  the 
mechanism  for  the  generation  of  turbulence  is  similar  in  cylindrical  and  planar  boundary  layers.31  This 
mechanism  is  the  burst  cycle3^  where  low-speed  fluid  intermittently  lifts  up  Irom  the  wall  resulting  in  an 
unstable  shear  layer  that  violently  breaks  up  generating  turbulent  velocity  fluctuations  followed  by  a 
sweep  of  fluid  toward  the  wail. 

Willmarth  et  al.12  recall  a  proposal  for  a  mechanism  for  the  cyclic  occurence  of  bursts  in  a 
planar  boundary  layer  that  suggests  that  the  pressure  fields  resulting  from  large  eddies  passing  over 
the  sublayer  produce  a  massaging  action  in  a  localized  region  that  intermittently  results  in  a  burst.  Since 
the  convection  velocity  in  a  cylindrical  boundary  iayer  is  ihe  same  as  for  a  planar  boundary  in  spite  of  the 
"fullness"  of  the  cylindrical  profile,  the  time  scale  of  the  massaging  action  is  similar  in  a  cylindrical 
boundary  layer.  Lueptow  and  Haritonidis31  support  the  existence  of  this  massaging  action  as 
evidenced  by  flow  visualization  of  large  scale  coherent  structures  and  velocity  correlations  between 
probes  near  the  wall  and  further  out  in  the  boundary  layer.  They  suggest  that  the  outer  fluid  "may  wash 
over  the  cylinder  precipitating  events  on  the  sides  of  the  cylinder."  Luxton  et  al.9  note  that  these 
washes  of  outer  fluid  may  sweep  away  the  inner  layer  maintaining  vorticity  in  the  layer  very  efficiently. 
Since  the  mixing  in  a  cylindrical  boundary  layer  is  controlled  by  the  inertia  of  the  large  eddies  rather  than 
by  viscous  effects  associated  with  the  wall,  the  profile  of  the  mean  velocity,  Reynolds  stress,  and 
intemnittency  are  altered  from  those  of  a  planar  boundary  layer. 
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6.  FURTHER  RESEARCH 


As  evident  throughout  this  report,  the  character  of  cylindrical  boundary  layers  is  poorly 
understood  in  comparison  to  planar  boundary  layers.  Several  potentially  productive  areas  of  research 
ara  outlined  below. 


*  MEASUREMENT  OF  MEAI I  VELOCITY  PROFILES 

-  Mean  velocity  data  are  needed  over  a  wider  range  of  conditions. 

In  spite  of  the  large  number  of  measurements  of  the  mean  velocity  profile,  a  definitive  scaling 
law  is  elusive.  This  is  a  result  of  inference  of  scaling  laws  from  poor  experimental  data  and  a  lack  of 
adequate  theoretical  models.  For  the  most  part,  though,  more  experimental  data  over  a  wide  range  of 
experimental  conditions  is  necessary.  For  instance,  the  data  of  Willmarth  et  at.12  appears  to  cover  a 
very  wide  range  of  nondimensional  parameters.  But  thL  variation  in  parameters  was  achieved  by  the 
variation  of  the  cylinder  diameter  and,  to  some  extent,  the  free  stream  velocity.  The  streamwise  location 
of  the  measurement  for  a  given  cylinder  was  not  an  independent  variable.  Ideally,  mean  velocity  profile 
measurements  should  be  made  while  varying  the  cylinder  diameter  (a).  Ihe  free  stream  velocity  (U„J.  and 
the  streamwise  measurement  location  (x).  As  Willmarth  et  al.1 2  point  out,  ‘...experiments  shou'd  be 
designed  to  provide  sufficient  data  concerning  the  effects  of  transverse  curvature  to  allow  a  complete 
empirical  formulation  of  the  mean  flow  field.  This  was  possible  in  the  tr/o-dimensiona!  [planar]  only  after 
sufficient  data  and  understanding  had  accumulated."  Further  analysis  of  the  presently  available  data 
may  also  be  helpful. 


WALL  SHEAR  STRESS 


The  skin  friction  in  a  cylindrical  boundary  layer  continues  to  be  an  unresolved  issue.  Sinca 
Richmond’s  work7  no  attempts  at  the  direct  measurement  of  the  coefficient  of  friction  have  been  made. 

It  seems  likely  that  tecnnological  advances  since  then  may  permit  such  measurements. 

-  Compilation  of  drag  measurement;, 

Many  measurements  have  been  mado  on  the  drag  of  cylinders  and  fibers.  These  results  could 
be  compi'ed  in  a  useful  format  to  provide  insignt  for  an  v— .pirical  relation  for  the  skin  friction  on  a  cylinder 
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-  Fluctuating  wall  shear  stress 

No  measurements  of  the  fluctuating  wall  shear  stress  on  a  cylinder  have  been  made. 
Measurements  using  hot  film  shear  stress  probes  should  be  made  on  cylinders  of  different  diameters 
under  varying  flow  conditions  to  fill  this  void. 

-  Drag  reduction 

Techniques  for  drag  reduction,  such  as  large  eddy  break-up  devices,  streamwise  riblets,  or 
microbubble  injection,  could  be  tested  in  a  cylindrical  boundary  layer.  The  cylindrical  geometry  has  the 
added  advantage  of  eliminating  edge  effects  that  are  present  in  a  planar  boundary  layer  experiment. 

•  MEASUREMENT  of  turbulent  quantities 

-  Measurements  of  the  Reynolds  stress,  turbulence  intensity,  and  velocily  spectra 

Measurements  of  the  turbulent  quantities  in  the  boundary  layer  are  more  scarce  than 
measurements  of  the  mean  velocity  profile.  Aside  from  measurement  difficulties,  the  streamwise  and 
wan-normal  turbulence  intensities  as  weii  as  the  Reynolds  stress  have  been  measured  over  oniy  a  very 
limited  range  of  parameters.  Extending  those  measurements  over  a  wider  range  of  parameters  should 
provide  needed  insight  into  the  character  of  a  cylindrical  boundary  layer 

-  Wall  pressure  fluctuations  for  8/a>4 

Wall  pressure  fluctuations  have  only  bean  measured  for  relatively  large  cylinders.  Using 
miniature  hearing  aid  microphones  the  wall  pressure  could  be  measured  for  much  smaller  cylinders. 
Arrays  of  microphones  could  be  used  to  map  contours  of  constant  wall  pressure  correlations  and  to 
measure  the  low  wavenumbor  portion  of  the  wall  pressure  spectrum. 

•  FLOW  NOISE  REDUCTION 

--  Techniques  for  reduction  of  low  frequency  velocity  fluctuations 

Lehman41  noted  a  reduction  in  the  kinetic  energy  of  the  velocity  fluctuations  at  low  wave- 
numbers  when  the  boundary  layer  was  exposed  to  a  local  transverse  motion  (S/a=0.2).  This  suggests 
that  further  investigation  of  a  method  tor  the  reduction  of  wall  pressure  fluctuations  based  on  local 
transverse  motion  of  the  wall  may  be  worthwhile. 
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•  STRUCTURE  OF  TURBULENCE 

-  Flow  visualization  to  character  outer  structures 

Flow  visualization  of  cylinders  of  different  diameters  moving  through  a  tank  of  quiescent  water 
may  provide  the  details  of  the  burst  cycle  and  the  effect  of  outer  eddies  on  the  generation  of  turbulence 
at  the  wall  for  cases  where  the  boundary  layer  thickness  is  much  greater  than  the  cylinder  radius. 

-  Space-time  correlations  of  velocity  to  characterize  coherent  structures 

Measurements  of  the  space-time  correlations  of  velocity  and  wall  pressure  for  cylinders  of 
different  diameters  may  provide  an  understanding  of  coherent  structures  in  the  flow  field.  Further  event 
detection  measurements  including  the*  measurement  of  the  correlation  between  wall  pressure  events 
and  velocity  events  may  be  helpful  in  providing  insight  into  the  wall  pressure  spectrum  and  its 
relationship  to  the  burst  cycle. 

•  MODELING  OF  CYLINDRIC/L  BOUNDARY  LAYERS 

•  .  -  Enhanced  models  of  the  integral  thicknesses  arid  skin  friction 

Past  predictive  models  of  the  coefficient  of  friction  have  been  band  on  very  crude  models  of 
the  boundary  layer  profile.  These  models  could  be  improved  and  enhanced  based  on  recent 
experimental  data.  Extension  of  these  techniques  to  rough  surfaces  may  be  useful. 


-  Direct  computer  simulation  of  a  cylindrical  boundary  layer 


Advances  in  computational  fluid  dynamics  may  permit  the  direct  computational  simulation  of  a 
turbulent  boundary  layer  on  a  cylinder.  The  turbulent  boundary  layer  ori  a  flat  plate  has  been 
successfully  modeled  at  a  Reynolds  number  based  on  momentum  thickness  of  up  to  14io.8^  Similar 
techniques  could  be  used  for  a  cylindrical  boundary  layer.  Commercially  available  computational  fluid 
dyriamics  codes  run  on  a  supercomputer  may  provide  useful  first-order  appioximations  to  the  flow. 

•  NON-AXISYMMETRIC  BODIES 

--  Effect  of  oval  cross-section  on  a  turbulent  boundary  layer 
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While  most  of  the  work  described  above  has  emphasized  the  boundary  layer  that  develops  on 
a  cylinder,  transverse  curvature  that  is  not  axisymmetric  has  not  been  studied.  For  instance,  the  radius 
of  transverse  curvatu-e  at  the  ends  of  the  major  axis  of  a  long,  oval  cross-section  body  will  be  much  less 
than  that  at  the  ends  of  the  minor  axis.  This  may  allow  the  isolation  of  the  effects  of  the  outer  flow  on  the 
generation  of  turbulence  from  the  effects  of  transverse  curvature,  since  large-scale  stmctures  cannot 
pass  as  easily  from  one  side  of  the  oval  cross-section  body  to  the  other  side. 
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